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PREFACE 


Tills  Memorandum  is  the  result  of  a  study  carried  out  at  RAHD  dur¬ 
ing  the  summer  of  1965-  The  objective  was  two-fold:  to  assess  the  ef¬ 
ficacy  of  the  theory  existing  at  that  time,  and  to  brief  RAND  personnel 
on  a  branch  of  fluid  mechanics  which  is  ordinarily  overlooked  in  texts 
and  graduate  curricula.  A  short  series  of  lectures  accompanied  the  de¬ 
velopment  of  this  Memorandum.  The  response  to  this  indicated  that  pub¬ 
lication  of  this  Memorandum  would  be  worthwhile  as  part  of  a  continuing 
RAND  program,  under  U.S.,  Air  Force  Project  RAND,  of  surveying  the  state 
of  the  art  in  aerodynamics. 

A  framework  for  the  description  of  three-dimensional  flow  fields 
is  a  requisite  for  the  rational  analysis  and  design  of  aerodynamic  sys¬ 
tems.  This  Memorandum  focuses  on  one  topic  in  this  general  area,  three- 
dimensional  boundary  layer  phenomena,  and  introduces  the  concepts  and 
techniques  of  this  subject  to  those  workers  in  aerodynamics  and  fluid 
mechanics  who  are  not  yet  actively  engaged  in  its  study. 

The  author  is  a  professor  in  the  Division  of  Aeronautical  Sciences 
at  the  University  of  California,  Berkeley,  and  a  consultant  to  The  RAND 
Corporation. 
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SUHMARY 

The  basic  concept  of  a  three-dimensional  boundary  layer  is  re¬ 
viewed  from  both  physical  and  mathematical  viewpoints-  Particular  em¬ 
phasis  is  placed  on  the  various  causes  of  secondary  flow,  with  geodesic 
curvature  of  the  surface  streamlines  of  inviscid  flow  receiving  the 
most  detailed  consideration.  Various  exact  solutions  for  steady,  in¬ 
compressible  laminar  flow  are  reviewed  and  a  proposal  for  a  finite- 
difference  scheme  for  arbitrary  inviscid  flows  and  surface  conditions 
is  sketched  (but  not  developed  and  tested).  The  momentum-integral 
method  and  other  approximation  schemes  are  briefly  discussed • 

Compressibility  effects  are  discussed  qualitatively,  and  then  at¬ 
tention  is  turned  to  the  stability  of  laminar  flows  and  the  transition 
to  turbulence. 

The  equations  for  time-averaged  turbulent  flow  are  derived  and 
criticised,  both  in  differential  and  integral  form.  We  review  the  two- 
dimensional  case  as  a  reminder  of  the  essential  difficulty  of  indeter¬ 
minacy,  and  then  critically  examine  existing  empirically-based  models 
of  the  three-dimensional  turbulent  boundary  layer.  No  attempt  is  made 
to  review  existing  methods  of  computing  the  development  of  turbulent 
layers . 

We  conclude  with  suggestions  for  further  study,  most  of  which  seem 
as  interesting  today  as  they  did  in  1965* 


S 
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SYMBOLS 


All  symbols  are  defined  where  they  appear  in  the  text,  and  many 
appear  in  only  one  section  of  the  Memorandum.  Some  symbols  are  used 
with  more  than  one  meaning  when  their  appearances  are  well  separated 
in  space  and  context,  and  confusion  seems  unlikely. 

The  symbols  which  recur  throughout  the  text  are  listed  below. 


Independent  Variables  and  Given  Parameters 


S,  "H 
C 

hl,h2 

H1 


2 


R 

L 

U,  V 

P 

A 

P 

P 

P 

v 


x,y,z 


ra,n,r,B 


coordinates  of  an  orthogonal  mesh  covering  the  body.  De¬ 
fined  in  dimensions  of  length 
distance  normal  to  the  body 

metric  coefficients  for  §  and  1)  (dimensionless) 

geodesic  curvature  of  the  coordinate  curve,  f]  =  constant, 
on  the  body  (length 

geodesic  curvature  of  §  -  constant  curve  (length  *) 

angular  velocity  of  reference  frame  relative  to  inertial 
space;  u>2,  are  its  §-,  T|-,  and  ^-components 

distance  from  the  point  (§,!], C)  t0  the  axis  of  reference 
frame  rotation 

a  characteristic  body  length 

I-  and  T|-components  of  velocity  at  the  wall  in  ir.viscid  flow 
pressure,  predetermined  by  the  inviscid  flow 
potential  of  a  conservative  body  force 
density,  assumed  constant  except  in  Section  XI 

2  2 

abbreviation  for  p  +  p(A  -  to  R  /2) 
viscosity,  assumed  constant  except  in  Section  XI 
kinematic  viscosity,  assumed  constant  except  in  Section  XI 
dimensionless,  and  in  the  cast  of  z,  stretched,  versions  of 
S,  T),  Q,  defined  on  page  13 

dimensionless  pressure-gradient  parameters  defined  on  pages 
13  and  14 


H,  X 


=  dimensioniess  coordinate-curvature  parameters  (proportional 
to  respectively),  defined  on  page  14 


Dependent  Variables 

•i,v,w  -  §- ,  tu  ,  ^-components  of  velocity  in  the  boundary  layer 

t'.g'  =  dimensionless  velocity  components  f'  =  u/U,  g'  =  v/V 

G^  =  dimensionless  v-velocity  in  intrinsic  coordinate  system, 
G '  =  u/U 

Q  -  vorticity  vector,  with  §- ,  "H- ,  and  ^-components 

=  surface  shear  stress  (skin  friction  vector),  with  §-  and 


"'i-components  r^, 


I.  INTRODUCTION 


When  one  considers  the  high  degree  of  three-dimensionality  char¬ 
acteristic  of  almost  all  bodies  of  aerodynamic  interest,  it  seems  re¬ 
markable  that  the  study  of  three-dimensional  boundary  layers  has  gained 
so  little  popularity  among  aerodynamicists.  Physico-chemical  or  elec¬ 
tromagnetic  complications  are  frequently  introduced  into  studies  of  two- 
dimensional  or  axisymmetric  boundary  layers,  but  geometrical  complica¬ 
tions  are  regarded  with  great  trepidation — if  regarded  at  all. 

The  possibility  that  this  attitude  not  only  defied  practical  real¬ 
ity  but  also  hampered  our  efforts  to  understand  fluid  motions  was  ini¬ 
tially  brought  to  my  attention  by  a  lecture  given  by  Dr.  E.  A.  Eichel- 
brenner  of  O.N.E.R.A. ,  by  Professor  M.  J.  Lighthill's  stimulating  chap¬ 
ter  in  Rosen'nead's  Laminar  Boundary  Layers,  and  by  a  great  number  of 
articles  in  the  meteorological  literature,  especially  those  describing 
attempts  to  integrate  the  "primitive  equation"  as  a  model  cf  large- 
scale  atmospheric  circulations. 

The  opportunity  to  investigate  the  matter  further  was  provided  by 
RAND,  at  the  suggestion  of  Dr.  Jerome  Aroesty.  Ky  early  impressions, 
gained  from  reading  the  survey  articles  listed  in  the  bibliography, 
have  confirmed  my  feeling  that  this  topic  is  worthy  of  a  more  broadly 
distributed  interest  among  fluid  dynamicists. 

The  objective  of  this  Memorandum  is  to  provide  a  nonspecialist's 
view  of  the  work  of  specialists.  It  is  assumed  that  the  reader  is  sub¬ 
stantially  familiar  with  two-dimensional  boundary  layer  theory,  and 
thus  has  a  feeling  for  how  the  magnitude  of  the  velocity  vector  is  re¬ 
duced  through  the  boundary  layer  under  various  sorts  of  external  flows. 
Our  discussions  of  three-dimensional  boundary  layers  will  focus  on  how 
the  direction  of  the  velocity  vector  varies  through  the  boundary  layer, 
and  hence  little  attention  will  be  paid  to  those  flows  (i.e. ,  axisym¬ 
metric,  nonspinning)  in  which  no  such  variation  appears. 

It  is  in  this  turning  of  the  velocity  vector  that  we  discover  the 
phenomenological  richness  of  three-dimensional  boundary- layer  flows. 
Freed  from  the  constraint  of  two-dimensional  motion,  the  fluid  seems 
to  "come  alive.,"  seeking  more  or  less  tangential  detours  around 
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"obstacles"  imposed  by  the  pressure  distribution,  bringing  about  a  fas¬ 
cinating  topology  of  skin  friction  lines  or  "surface  streamlines,"  com¬ 
plicating  the  concept  cf  separation,  and  so  on.  We  would  hope  to 
achieve  from  our  study  some  feeling  for  how  and  why  these  things  occur, 
perhaps  some  rules  of  thumb  concerning  the  effects  of  body  spin  or  of 
inviscid  streamline  curvature  to  add  to  those  for  the  effects  of  stream- 
wise  pressure  gradient. 

It  should  be  understood  at  the  outset  that  this  Memorandum  does 
not  constitute  a  review  article  in  the  usual  sense-  In  the  first 
place  T.  am  not  personally  expert  in  the  subject  treated,  and  have  de¬ 
rived  my  information  entirely  from  the  writing  of  others-  Topics  of 
real  practical  importance,  such  as  the  effects  of  compressibility  and 
of  spinning  bcdies,  have  frankly  not  been  accorded  their  deserved 
share  of  attention,  partly  by  intent  (to  focus  attention  more  sharply 
on  purely  geometrical  complexities)  and  partly  by  the  limitations  of 
time  and  my  acquaintance  with  the  literature. 

This  Memorandum  begins  with  a  largely  intuitive  review  (and  exten¬ 
sion  to  three-dimensional  flow)  of  the  concept  of  a  boundary  layer  and 
of  a  properly  posed  problem  of  steady  flow  in  boundary- layer  theory. 

The  laminar  boundary- layer  equations  are  then  presented  in  orthogonal 
curvilinear  coordinates.  For  steady,  incompressible  flow  we  introduce 
the  Falkner  transformation  to  eliminate  the  normal  velocity  and  to 
state  the  tangential  momentum  equations  in  a  form  convenient  f*„r  all 
our  future  discussions  of  laminar  flow.  These  discussions  begin  with 
a  suggestion  for  finite-difference  solution  procedures  in  various  co¬ 
ordinate  systems,  following  the  method  applied  to  two-dimensional  flows 
by  Smith  and  Clutter  (1963),  but  making  frequent  reference  to  the 
method  developed  by  Raetz  (1957)  for  three-dimensional  flows.  Efforts 
to  depart  from  the  same  starting  point  by  series  expansion  in  the  sur¬ 
face  coordinates  are  noted.  A  brief  review  of  exact  similar  solutions 
at  stagnation  points  and  for  flow  over  developable  surfaces  follows. 

Next  treated  are  the  two  principles  that  have  been  exploited  in 
most  solutions  found  to  date,  namely  the  Independence  principle  for 
swept  flows  and  the  prevalence  principle  for  cases  with  weak  secondary 
flow. 
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Returning  to  a  more  general  viewpoint,  we  give  a  semi-intuitive 
discussion  of  the  causes  and  nature  of  secondary  flows,  and  a  fairly 
extended  discussion  of  the  topology  of  skin  friction  lines,  or  "sur- 
face  streamlines,”  concluding  with  a  definition  and  discussion  of 
boundary- layer  separation  in  three  dimensions. 

Momentum-integral  approximation  methods  are  discussed  next,  both 
because  they  lead  to  some  useful  results  with  a  considerable  savings 
in  computational  effort,  and  because  the  difficulties  into  which  they 
occasionally  lead  us  provide  useful  tests  cf  our  intuitive  understand¬ 
ing  of  the  boundary  layer.  The  momentum- integral  equations  are  de¬ 
rived  directly  from  the  Fa lkner- transformed  equations  for  laminar  flew, 
and  are  later  re-derived  more  conventionally  for  turbulent  flow,  to 
provide  a  check.  Velocity  profile  assumptions  and  the  use  of  weighted 
integral  equations  are  discussed  at  some  length. 

Compressibility  effects  are  mentioned  very  briefly  and  mostly  by 
reference  in  Section  XI,  as  are  laminar-flow  stability  in  Section  XII 
and  transition  to  turbulence  in  Section  XIII.  The  displacement  thick¬ 
ness  of  a  three-dimensional  boundary  layer  is  defined  in  Section  XI. 

An  extended  discussion  of  the  incompressible  turbulent  boundary 
layer  follows  (Section  XIV),  with  the  assumption  that  the  reader  can 
benefit  from  a  brief  review  of  the  problem  of  formulating  a  momentum- 
integral  approach  for  two-dimensional  mean  flows.  It  is  concluded 
that  very  successful  three-parameter  models  of  mean  velocity  profiles 
exist,  but  that  methods  for  predicting  the  variation  of  the  profile 
parameters  are  still  essentially  deficient.  No  three-parameter  pro¬ 
file  model  as  generally  successful  as  the  '‘wall-wake”  model  of  two- 
dimensional  flow  has  yet  been  found  for  the  three-dimensional  case. 

Finally,  the  Memorandum  concludes  with  a  number  of  suggestions 
for  experiments,  calculations,  or  analyses  which  appear  potentially 
fruitful  for  developing  cur  powers  of  description  and  prediction  in 
the  areas  covered. 

The  bibliograpy  is  not  complete,  but  does  include  all  available 
recent  survey  articles.  These  in  turn  contain  comprehensive  biblio¬ 
graphies,  from  which  almost  all  the  individual  references  given  here 
are  drawn.  One  gets  the  impression  from  current  issues  of  abstracting 


journals  that  the  subject  currently  receives  little  attention  in  com¬ 
parison  with  more  recently  glamorous  subjects,  so  that  even  a  1962  or 
1963  bibliography  or  survey  article  may  be  considered  reasonably  up 
to  date.  The  sources  which  have  most  influenced  the  spirit  of  this 
Memorandum  are  the  survey  articles  by  Mager  (1964),  Head  (1961),  and 
Rotta  (1962),  the  chapters  by  Lighthill  (1963),  the  report  by  Raetz 
(1957),  the  paper  of  Smith  and  Clutter  (1963),  the  N. L.L.  reports  of 
Tinman  (1951)  and  Zaat  (1956),  the  paper  of  Coles  (1956)  and  the 
O.N.E.R.A.  publications  and  other  papers  by  Eichelbrenner  and  co¬ 
workers  . 
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II.  THE  BOUNDARY - LAYER  APPROXIMATION  IN  THREE-DIMENSIONAL  FLOW 


Survey  articles  (e.g.,  Moore  (1956)  or  Mager  (1963))  usually  dis¬ 
tinguish  between  boundary- layer  flows  and  boundary-region  flows.  In 
the  examples  usually  cited,  the  distinction  seems  to  be  the  following: 

(A) .  In  boundary- layer  flows  we  may  treat  diffusion  as  pro¬ 

ceeding  in  only  one  direction,  normal  to  the  wall. 

This  generalizes  the  concepts  of  two-dimensional  boundary- layer 
theory,  in  which  streanrwise  diffusion  is  neglected,  to  state  that  dif¬ 
fusion  in  all  directions  parallel  to  the  wall  is  to  be  neglected. 

The  physical  idea  behind  this  approximation  is  the  thinness  of 

the  boundary  layer,  as  measured  along  a  normal  to  the  wall,  compared 

* 

to  typical  streamwise  or  lateral  distances  of  interest. 

In  boundary-region  flows,  exemplified  by  flow  along  a  corner  or 
edge  where  two  walls  join  abruptly,  there  is  clearly  no  unique  normal 
direction  at  a  boundary- layer  point  close  to  the  corner,  and  hence 
lateral  diffusion  must  be  reinstated.  These  flows  are  not  treated 
here- 

(B) .  In  boundary- layer  flows  we  may  take  the  pressure  to  be 

predetermined  (on  first  approximation)  by  the  invisci'.- 
flow  pattern  on  the  s'^face  cf  the  given  pocy. 

The  accuracy  of  this  assumption  depends  upon  (I)  the  boundary 
layer  being  thin  relative  to  any  local  radius  of  curvature  of  the  body, 
and  (2)  the  "displacement  effect"  of  the  boundary  ?.ayer  producing  only 
small  perturbations  in  the  local  external  flow. 

Both  assumptions  A  and  B  may  be  violated  by  th'  phenomenon  of  sep¬ 
aration,  if  this  leads  to  detachment  of  a  thin  shear  layer  the  normal 
to  which  does  not  nearly  coincide  with  the  local  norma-  to  the  wall. 

In  such  a  case  assumption  A  would  misdirect  the  dominant  diffusive 
flux,  and  the  large  displacement  effect  would  violate  assumption  B. 


We  shall  often  need  a  word  to  indicate  the  direction,  parallel 
to  the  wall,  but  normal  to  the  local  velocity  vector.  For  this  we 
choose  "lateral." 
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III t  BOPKDARf  VALUES  AND  INITIAL  DATA 

tne  specification  of  boundary  valueB  for  three-dimensional  bound¬ 
ary-layer  problems  exactly  parallels  that  for  two-dimensional  flows. 

We  ordinarily  specify  no  slip  at  the  wall  for  the  tangential  compo¬ 
nents  of  velocity,  and  may  prescribe  an  arbitrary  distribution  of  nor¬ 
mal  velocity  (suction  or  blowing)  on  the  wall,  so  long  as  we  do  not 
thereby  excessively  thicken  the  boundary  layer.  (Roughly  speaking, 
these  normal  velocities  should  only  be  of  the  order  of  magnitude  of 
those  arising  spontaneously  in  the  boundary  layer  in  the  impermeable 
wall  case.)  Temperature  or  heat  transfer,  and  concentration  or  mass 
transfer  must  also  be  given  at  the  wall. 

At  "infinity11  with  respect  to  a  suitably  stretched  normal  distance 
variable,  the  tangential  velocity  components ,  temperature  and  composi¬ 
tion  arc  '«p«ated  to  the  re Lues  of  these  quantities  appearing  at  tha 
wall  in  the  inviscid  flow  solution.  The  normal  velocity  is  not  speci¬ 
fied  there.,  but  is  derived  from  tha  boundary- layer  solution. 

As  v?  shall  discuss  in  more  detail  in  Section  VIII,  the  boundary 
layer  in  steady  flow  over  a  smoothly  curved  body  can  be  said  to  be 
born  at  a  finite  number  of  isolated  nodal  points  of  attachment,  some 
of  which  coincide  with  the  attachment  points,  or  "forward  stagnation 
points,"  of  the  inviscid  flow.  In  the  immediate  vicinity  of  such  at¬ 
tachment  points  can  be  found  a  local  similar  solution  of  the  boundary- 
layer  equations,  which  is  also  an  exact  solution  of  the  full  Navier- 
Stokes  equations  in  an  infinitesimal  neighborhood  of  the  normal  at  the 
stagnation  point.  These  special  solutions  can  be  used  to  establish 
initial  data  for  a  boundary- layer  calculation  in  an  extended  region 
surrounding  the  attachment  point  by  assigning  values  to  all  unknowns 
on  a  "wall"  of  normals  arising  out  of  a  conveniently  chosen  surface 
curve  enclosing  the  attachment  point. 

If  the  inviscid  flow  attachment  is  along  an  infinitely  sharp  lead¬ 
ing  edge  (or  an  edge  which  is  conveniently  idealised  as  such)  then  the 
boundary  layer  is  born  all  along  this  edge  and  its  initial  development 
is  given  by  a  similar  solution  «.f  a  Falkner-Skan  type.  This  can  be 
used  to  construct  initial  data  on  a  normal  wall  closely  paralleling 
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the  leading  edge.  This  procedure  would  also  apply  to  the  degenerate 
case  of  essentially  two-dimensional  flow  in  which  the  leading  edge  is 
a  cylinder  of  finite  curvature  aligned  exactly  normal  to  the  oncoming 
flow.  (Think,  for  example,  of  a  ring  wing  or  axisynsnetric  engine  in¬ 
take  at  zero  incidence.) 

So  far,  the  discussion  resembles  the  corresponding  one  of  two- 
dimensional  boundary  layers.  Complications  enter,  however,  when  we 
consider  the  effects  of  three-dimensionality  on  the  domain  of  depen¬ 
dence  downstream  of  a  finite  segment  of  our  "initial  data  curve." 

Having  adopted  the  assumption  that,  with  boundary  conditions 
given  as  above,  the  boundary  layer  evolves  under  the  influence  of  a 
predetermined  pressure  field  in  a  region  in  which  diffusion  takes  place 
only  along  the  normal  to  the  wall,  we  expect  to  be  able  to  construct 
solutions  by  starting  with  the  given  values  on  the  "initial  data  sur¬ 
face"  and  "marching"  sccadily  downstream,  accounting  at  each  step  for 
the  modification  of  the  velocities,  etc.,  by  convection  along  the  lo¬ 
cal  velocity  vector  and  diffusion  along  the  local  wall  normal.  An  im¬ 
mediately  practical  question  is  the  following:  "Over  what  area  of  the 
body  surface  will  the  solution  be  completely  determined  by  the  given 
pressure  distribution  and  inviscid  velocity  field,  and  by  the  solution 
already  known  over  a  given  finite  segment  of  the  'initial  data  curve'?" 

In  Fig.  1,  AB  is  the  finite  segment  in  question.  We  assume  that 
Ab,  and  the  entire  curve  of  which  it  is  part,  is  oriented  so  that  all 
boundary- layer  particles  cross  it  from  one  side  to  the  other.  This 
condition,  which  physically  requires  that  all  flow  entering  the  region 
downstream  of  the  "initial  data  surface"  comes  from  a  single  "boundary- 
layer  birthplace"  or  attachment  region,  amounts  mathematically  to  the 
requirement  that  the  initial  data  surface  not  be  tangent  to,  or  crossed 
more  than  once  by,  a  characteristic  curve  of  the  boundary- layer  equa- 
tione  (see  Raetz  (1957)  for  mathematical  discussions). 

* 

Raetz  states  an  "influence  principle"  as  follows.  "The  influ¬ 
ence  of  the  solution  at  any  point  is  transferred  to  other  points  first 
by  conduction  along  the  straight  line  paralleling  the  £  axis  [the  nor¬ 
mal]  and  passing  through  that  point  and  then  by  convection  downstream 
along  all  streamlines  through  that  line."  He  calls  the  dashed  lines 
through  A  and  B  the  "outer  and  inner  characteristic  envelopes." 
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Fig.  1  —  Regions  of  influence  and  domain  of  dependence. 


Particles  which  follow  varied  downstream  trajectories  will  pass 
through  any  of  the  normals  generating  the  surface  over  AB,  according 
to  their  initial  distance  from  the  wall. 

If  we  imagine  the  normals  through  A  and  B  to  be  line  sources  of 
dye,  the  plan  view  of  the  distorting  dye  sheets  in  the  downstream  re¬ 
gion  would  exhibit  the  shaded  banners  of  the  sketch. 

According  to  our  assumption  of  purely  normal  diffusion,  every 
particle  in  the  boundary  layer  over  these  shaded  banners  can  become 
stained  by  this  dye,  whether  or  not  it  passed  through  the  normal  over 
A  or  B,  but  no  other  particles  are  contaminated.  The  shaded  regions 
can  then  be  called  the  regions  of  influence  of  the  initial  data  lines 
through  A  and  B. 

Correspondingly,  the  boundary  layer  over  the  unshaded  region  be¬ 
tween  the  two  banners  knows  only  about  the  past  history  of  particles 
which  have  passed  the  initial  data  surface  over  the  arc  AB,  and  hence 
the  unshaded  area  is  the  area  referred  to  in  the  underlined  question 
above.  We  shall  call  it  the  domain  of  dependence  of  the  solution  over 
the  arc  AB,  although  in  the  strict  mathematical  sense  that  term  should 
include  also  the  region  upstream  of  AB,  defined  by  stream  surfaces 
which  eventually  converge  at  the  attachment  normal. 
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It  may  appear  startling  to  assert  that  once  we  have  adopted  the 
boundary- layer  approximations  we  have,  in  effect,  no  further  need  of 
the  fluid  which  passes  outboard  of  the  lines  over  A  and  B,  in  order 
to  determine  the  flow  in  our  unshaded  region.1  However,  there  are  but 
two  ways  in  which  this  fluid  can  make  its  presence  felt.  The  first  is 
by  pressure,  but  that  should  be  predetermined.  The  second  is  by  lat¬ 
eral  diffusion,  but  we  have  ignored  that  for  the  moment. 
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IV.  BOUNDARY- LAYER  EQUATIONS  IK  ORTHOGONAL  COORDINATES :  LAMINAR  FLCM 

CONTINUITY  AND  MOMENTUM  EQUATIONS 

Following  the  notation  of  Mager  (1963),  to  whom  the  reader  is  re¬ 
ferred  for  detailed  discussion  and  references,  we  introduce  general¬ 
ized  orthogonal  coordinates  §,  "H,  £.  The  lines  of  constant  §  and  T) 
form  a  network  on  the  body  surface,  and  £  increases  away  from  the  body 
along  the  normal. 

The  increment  of  distance  between  coordinate  lines  §  and  §  +  d§ 
is  hj  df,;  corresponding  increments  are  df]  and  h^  d£.  For  thin 
boundary  layers  we  can  assume  that  the  metric  coefficients  satisfy 

hl  =  h2  =  “  1 

so  that  £  can  be  interpreted  as  actual  normal  distance  from  the  wall, 

* 

which  lies  at  £  =>  0. 

The  velocity  components  in  the  §,  'll,  and  £  directions  are  called 
u,  v,  and  w.  The  corresponding  surface  velocities  of  the  inviscid 
flow  are 


U(5,T1),  V(§,Tp,  and  0 

The  boundary- layer  assumptions  reduce  the  complex  expressions  for 
viscous  stress  components  in  a  general  orthogonal  coordinate  system  to 
the  two  simple  results 

,  „  <3u  &v 

T31  ■  ac*  T32  “  ^  ac 


We  do  not,  of  course,  have  to  define  £  as  the  normal  distance 
itself.  Raetz,  for  example,  takes 


£  -  (1  -  u/U)l/2 

after  demanding  that  §  and  T)  be  so  oriented  as  to  render  u/U  monotonic 
in  normal  distance. 
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The  continuity  equation  is  simply 

If  +  flc  (hgUp)  +  ^  +  |j  Chih2“p>}  -  o 


The  ^-component  or  the  momentum  equation  is 


,  1  9A  u ) 

+  T~  Tc  ~  2uj-V  -  TT- 


uv 

“"1 

*»  V! 

V 

hlh2 

91) 

-  hlh2  a 

.  J_ 

9p 

1  9 

Ph'i 

IX 

P  dC 

The  11-component  reads 


5v  .  u  9v  .  v  9v  .  9v  u2  ^1  uv  ^‘2 

at  T  h,  9§  h,  OT]  w  dc  “  hxh2  afrj  h,h0  &§ 


12 


1_  9A 

h2  an 


,  .  tu2  a  d2  .  i  9p  i  a  z'  9v\ 

+  2(i)„u  -  -rr—  t— -  R  +  ~ —  xt 7  =*  —  r-r  \u  t— ) 

3  2h„  911  ph2  oil  P  9C  v  &C' 


The  ^-component  reads 

9A  ..  %  u)2  9  _2  1  Sp 

5c  '  2(  V  ‘  -  1  sc  R  +  5  ac  " 

In  the  §-  and  H-component  equation,  the  top  lines  simply  give  the 
components  of  acceleration,  while  the  terms  in  the  bottom  lines  repre¬ 
sent  (1)  a  body  force  derivable  from  the  potential  A,  (2)  a  Coriolis 
force  arising  if  the  body  (and  hence  the  coordinate  system,  is  spin¬ 
ning  uniformly  and  steadily  with  respect  to  an  inertial  frame,  (3)  a 
centrifugal  force  derived  from  such  spinning,  (4)  the  predetermined 
pressure  force,  and  (5)  the  viscous  force. 
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The  local  §,  "H,  and  Q  components  of  the  angular  velocity  of  the 

2  2  2  2 

body  are  -j^,  and  i^,  with  uT  =  +  X  +  u^*  The  distance  from 

the  spin  axis  to  the  point  (5,T|,Q  is  R. 

The  ^-component  equation  asserts  that  the  pressure  gradient  nor¬ 
mal  to  the  wall  counterbalances  the  body  force  and  the  "forces”  due 
to  spinning.  These  forces  are  assumed  to  be  sufficiently  small  so 
that  throughout  the  thin  boundary  layer  we  may  take 


?  =  P(§,T») 


with  p  relate*  to  the  velocities  U  and  V  by 


au  .  u  au  v  av  uv  Shi  v2  ah2 
at  •  hL  +  h2  an  +  hLh2  an  '  h^  a§ 


.  1  d/\  0  „ 

+  h[  al  ‘  2uj3v 


“L_  <L  Rz  +  _JL  *£ 
2h,  a§  K  p  h.  a§ 


_L  Ae  -  0 

ehl  « 


av  uav  vav  ahi  uv  5h2 

at  +  hxa§  h2an  "  h^  an  h^  ai 


+  i_  +  2us  u  -  — —  R2  +  — L  Ae  =  o 

h2  an  +  ZV  2h2  an  Peh2  an 


u)  a  n2 


Here,  designates  the  surface  density  of  the  inviscid  flow. 

In  compressible  flows,  in  which  p  and  p  must  be  treated  as  vari¬ 
able,  we  need  also  the  energy  equation,  an  equation  of  state,  and  a 
description  of  the  dependence  of  p  and  of  the  thermal  conductivity 
upon  the  thermodynamic  variables.  In  what  follows  immediately,  we 
shall  assume  p  =  p  =  constant,  and  p  =  constant. 

This  assumption  focuses  attention  on  effects  due  to  the  three- 
dimensionality  of  geometrical  constraints,  or  to  Coriolis  force.  When 
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p  -  o&,  gravity  and  centrifugal  forces  act  equally  on  boundary- layer 
and  external  flow,  and  hence  this  assumption  rules  out  qualitatively 
interesting  phenomena  such  as  free  convection.  We  shall  comment  brief¬ 
ly  on  these  and  other,  primarily  quantitative,  effects  of  variable  den¬ 
sity  and  viscosity  in  Section  XI. 

THE  FALKNER  TRANSFORMATION 

For  discussions  of  exact  self-similar  solutions  and  foi  the  for¬ 
mulation  of  systematic  methods  of  calculation  of  ,*.ore  genera^  boundary- 
layer  flows,  it  is  helpful  to  transform  variables  in  some  manner  which 
accounts  roughly  for  the  anticipated  magnitudes  of  ooundary- layer  thick¬ 
ness  and  velocity  components.  Many  such  transformations  are  available. 

We  choose,  rather  arbitrarily,  that  associated  with  the  name  of  Falkner,* 

We  call 


where  L  is  a  characteristic  body  dimension  and  v  is  the  kinematic  vis- 
cosity.  Thus 


_  Sy  _  1.  dx  _  dx  dy  _  dy 

~  an  l’  an  "  ac  ”  “  3C 


dz  _/  m  -  1\  z  dz  /n\  z_  dz  /~U~ 

\  2  I  xL’  dT|  \2/yL  *  dC  =  <*vLx 


where 


m(x,y; 


d  in  U 
d  in  x* 


n(x,y) 


d  in  U 
d  in  y 


* 

Laminar  Boundary  Layers,  L.  Rosenhead,  ed. ,  Oxford,  1963,  p. 

266  (hereafter  cited  in  the  text  as  L. B.L. ). 

irk 

The  choice  of  U  rather  than  V  and  of  §  rather  than  T\  in  the  nor¬ 
malization  of  Q  implies  only  that  §  increases  in  a  racre-or-less  down¬ 
stream  direction,  and  that  U  ^  0  except  at  singular  points  or  lines. 
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Later  we  shall  also  want  the  abbreviations 


/  x  _  o  in  V 

r(*’y>  -  jPZiTx’ 


s(x,y)  = 


5  /n  V 
S  in  y 


Our  dependent  variables  are  written  as 


u  =  U(x,y)  —  f(x,y,z) 


which  we  designate  simply  u  =  Uf '.  At  this  point  the  prime  of  f  does 
not  Imply  that  f  depends  only  upon  z,  but  is  simply  a  shorthand.  We 
shall  write  out  df/dx  and  df/dy  fully.  Similarly,  we  write 

v  h  V(x,y)  ~  g(x,y,z)  *  Vg' 

Using  the  continuity  equation,  with  p  »  constant,  we  find,  for  the 
simple  case  of  an  impervious  wall, 


+  vt(a  -f  +  08  +  fzs'  +  y§] } 


where  we  have  introduced  our  final  shorthand  notation 


9  In  h,  d  In  h^ 

H  dXnx’  ^  ~  5  in  y 


If  u,  v,  and  w  are  to  vanish  at  z  »  0,  we  see  that  boundary  con¬ 
ditions  on  f  and  g  are 


f(x,y,0)  >=*  f'(x,y,0)  «  g(x,y,0)  =•  g'(x,y,0)  -  0 
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while 


f'(x»y>®)  =  g'(x,y,®)  =  1 


The  transformed  component  momentum  equations,  in  which  the  pres¬ 
sure  gradients  have  been  eliminated  in  favor  of  inviscid-flow  velocity 
gradients,  etc-,  become 


rn,  .  1  /ra  +  1  ,  \  .  m  ,,  ^,2,  ,  Vx  1  17  n  .  , 

£  +ir  v- 2“  +  f£  +  ir(  >+^irUs  -rii8f 


and 


THE  TRANSFORMED  EQUATIONS  IN  "INTRINSIC  COORDINATES" 

Of  the  various  possible  special  orientations  ot  the  Q  and  1]  axes, 
one  of  the  most  interesting  is  that  in  which  §  measures  distance  down¬ 
stream  along  the  surface  streamlines  of  the  inviscid  flow.  The  corre- 

ie 

For  convenience,  all  numbered  equations  appearing  in  the  text 
are  collected  in  the  Appendix,  p.  89. 


-16- 


sponding  surface  coordinates  are  railed  "intrinsic  coordinates."  This 
choice  of  coordinate  system  implies  that  V  =  0,  and  calls  for  a  renor¬ 
malization  of  v,  as,  for  example, 

v  =  U(x,y)G  '(x,y,z) 

Then  G /(x,y,0)  =  G'CXjy,®)  =  0. 

The  transformed  momentum  equations  become  slightly  simpler,  namely 
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V.  EXACT  SOLUTIONS,  INDEPENDENCE  PRINCIPLE,  AND  CONDITIONS  FOR 
ABSENCE  OF  SECONDARY  FLOW 


STAGNATION  POINT  SOLUTIONS 

Suppose  we  deal  with  a  “rounded"  body,  to  which  the  flow  from  up¬ 
stream  "attaches"  at  some  point  P,  which  we  take  as  the  origin  of  co¬ 
ordinates.  In  the  immediate  vicinity  of  P  we  can  expand  the  surface 
metric  coefficients  in  power  series  in  x  and  y,  e.g.. 


hjU.y) 


h^O.O)  + 


(0,0)x  +  {^)  (0,0)y  +  ... 


In  particular,  we  choose  a  system  of  coordinates  that  is  locally  rec¬ 
tangular  at  the  origin,  so  that 

^(0,0)  =  h2(0,0)  =  I 


Since  (dh^/cbO^  ^  and  the  similar  derivatives  are  finite,  we  have, 
in  the  vicinity  of  the  origin. 


Furthermore,  if  the  approaching  upstream  flow  is  irrotational,  it 
can  be  shewn  that  there  is  an  orientation  of  the  x  and  y  axis  for  which, 
in  the  vicinity  of  the  origin, 

U  =*  Ax  +  . . . ,  V  *  By  +  . . . 

The  desired  orientation  of  the  x  and  y  axes  is  along  the  principal 
directions  of  curvature  of  the  surface.  To  be  specific,  we  may  take 
the  x  axis  along  the  direction  of  maximum  (convex  outward)  curvature. 
Then  we  shall  have  A  >  0  and  A  >  B. 

If  B  is  positive,  we  speak  of  a  nodal  attachment  point;  if  it  is 
negative,  we  designate  it  as  a  saddle  point  of  attachment. 
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We  see  chat  the  given  U  and  V  correspond  to 

m  =  8  =  1,  n  =  r  =  0 

All  the  parameters  on  the  left-hand  side  of  Eqs.  (1)  and  (2),  in¬ 
cluding  the  new  parameter 

_  _  B  Vx 

C  A  =  Uy 

are  constants.  We  exclude  spinning  of  the  coordinate  system,  so  oi^  =  0 

We  can  therefore  expect  to  find  similar  solutions  for  f  and  g,  for 
which  df/3x  =  dg/dx  =  0.  The  governing  equations  reduce  to 

f "'-i-  (f  +  Cg) f"  +  1  -  f /2  -  0 

g*'  +  (f  +  Cg)gv  +C(1  -  g/2)  -  0 

These  have  been  solved  numerically  by  Howarth  (1951)  for  0  s  C  i  1 
and  by  Davey  (1961)  for  -1  s  C  S  0.  Extension  to  compressible  flow  has 
been  made  by  Poots  (1965)  for  -.5  SCSI, 

The  limiting  cases  C  =  0  and  C  =  1  clearly  correspond  to  the  two- 
dimensional  and  axisymmetric  stagnation  points,  respectively. 

We  have  mentioned  the  practical  importance  of  the  nodal  solutions 
in  Section  III;  tabulation  and  further  discussion  are  given  on  pp. 

461  to  467  of  L. B.L.  The  saddle-point  solutions  are  particularly  in¬ 
teresting,  in  that  for  0  >  C  >  -0,4294  they  propose  a  locally  deter¬ 
mined  solution  for  a  region  of  flow  which  lies  within  the  domain  of 
dependence  of  initial  data  curves  surrounding  neighboring  nodal  points. 
This  is  not  permissible  in  principle,  but  may  be  approximately  useful 
in  practice  if  the  boundary  layer  has  effectively  "forgotten"  its 
initial  data  by  the  time  it  gets  close  to  the  saddle  point.  For 
-0.4294  >  C  >  -1,  Davey's  solutions  are  particularly  intriguing,  show¬ 
ing  a  streamline  pattern  which  changes  over  from  saddle-  to  node-like 
behavior  as  z  decreases  towards  zero.  The  g'  profile  shows  reverse 
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flow  near  the  wall  and  the  flow  thus  appears  in  a  sense  to  be  separated. 
(A  careful  discussion  of  separation  follows  in  a  later  section.)  Never¬ 
theless,  for  C  >  -1,  the  displacement  thickness  of  the  boundary  layer 
is  finite  and  the  external  flow  is  still  approaching  the  wall  along  the 
stagnation  streamline.  The  saddle-point  region  for  -0.4294  >  c  >  -1 
does  not  lie  within  the  domain  of  dependence  of  neighboring  nodal  points, 
and  in  many  ways  seems  to  be  a  region  wetted  not  directly  by  the  flow 
from  upstream,  but  by  a  very  simple  type  of  wake  flow. 

"DOUBLE"  FALKNER-SKAN  FLOWS  OVER  DEVELOPABLE  SURFACES 

Developable  surfaces  can  be  "unrolled"  into  a  plane  without  wrin¬ 
kling,  and  hence  our  (5,^)  coordinates  on  such  surfaces  can  be  every¬ 
where  rectangular.  This  makes 

hj  3  h£  *  1,  k  =  X  0  everywhere 

Now  if  **  0  (no  spinning)  and  if  we  can  find  an  orientation  -c  -  J 
y  axes  such  that 

it  »  m  n  ,  „  _  r  s 

U  *  Ax  y  and  V  «  Bx  y 

with  constant  A,  B,  m,  n,  r,  and  s,  and  if  the  constants  are  related 
by 

r  ■  m  -  1  and  s  *  n  +  1 

there  will  exist  similar  solutions  governed  by 

+  t?  +m(l  -  f'2)  gf'  +  n(l  -  g'f')}  «  0 

and 
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* 


These  have  been  solved  by  Yohner  and  Hansen  (1958)  for  B/A  =  1 
and  2,  and  for  many  combinations  of  m  and  n.  The  assumed  inviscid 
flows  are  in  general  quite  rotational,  and  the  resulting  nonuniformi- 
ties  of  total  pressure  lead  in  some  cases  to  "overshooting"  f  '-pro¬ 
files. 

Simi  lar  solutions  of  a  particularly  interesting  type  are  also  ob¬ 
tained  for  arbitrary  r  and  m,  providing  that  n  »  s  =  0.  The  inviscid 
flows  are  again  rotational  for  r  >  0,  and  the  governing  equations  are 

t"'  +  ff^  +  m(l  -  f  '^)  =  0  (The  Falkner-Skan  equation) 


and 


g"+  ®±I  fg'  +  r(1  .  f'g')  »  o 

These  have  been  solved  by  Hansen  and  Herzig  (1956)  for  m  =  0  and  for 
integer  values  of  r  from  1  to  10.  These  authors  noticed  that  the  equa¬ 
tion  for  z'  does  not  contain  g*,  and  that  the  equation  for  g r  is  linear, 
when  f '  is  considered  to  be  known.  Thus  solutions  for  any  fixed  m, 
but  with  V  given  by  a  polynomial  in  x,  can  be  generated  by  superposi¬ 
tion.  They  give  examples  from  the  field  of  turbomachinery  and  exhibit 
some  interesting  comparisons  with  experiment.  Yohner  and  Hansen  (1958) 
also  consider  this  case,  and  obtain  solutions  for  all  combinations  ob¬ 
tainable  from  m  =  0,  1,  2,  4,  6,  8,  10  and  r  =  0,  1,  2,  4,  6,  8,  10. 

THE  INDEPENDENCE  PRINCIPLE 

The  equations  treated  by  Hansen  and  Herzig  (1956)  exhibit  the  spe¬ 
cial  feature  that  f  can  be  found  independently  of  g,  and  we  may  inquire 
as  to  how  generally  this  can  be  done. 

We  examine  Eq.  (1)  and  observe  that  g  appears  only  in  terms  con¬ 
taining  derivatives  with  respect  to  y  (including  the  parameters  X,  n, 
and  s),  the  curvature  parameter  h,  and  the  coordinate  rotation  speed 
u )y  Thus,  in  a  y-independent  flow  over  a  nonspinning  body  with  a  de¬ 
velopable  surface,  f  is  independent  of  g.  This  is  the  independence 
principle. 
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A  sample  flow  would  be  that  over  a  yawed,  infinite  cylinder,  with 
y  measured  along  the  generators  of  the  cylinder.  Then  V  is  a  constant 
and  the  governing  equations  are 


The  equation  for  f  is  that  solved  by  Smith  and  Clutter  (1963a) 

for  a  variety  of  x-wise  pressure  distributions.  The  equation  for  g 

is  a  second-order  linear  equation  for  g so  that  the  generation  of 

solutions  for  the  spanwise  flow  (g  7)  for  any  of  Smith  and  Clutter's 

chordwise  flows  would  be  a  relatively  simple  matter.  A  procedure  for 

doing  this  numerically  is  suggested  by  Lindfield,  Pinsent,  and  Pin- 
* 

sent. 

For  the  special  case  m  =  constant,  the  equation  for  f  is  the 
Falkner-Skan  equation.  The  corresponding  similar  span-wise  profiles, 
g7,  have  been  found  by  J.  C.  Cooke  (1950),  whose  results  are  tabulated 
on  p-  471  of  L.B.L.  (In  that  table  they  are  called  g.)  Clearly  the 
case  m  =  1  gives  the  same  f  and  g  as  Howarth's  stagnation  point  solu¬ 
tion  for  C  =  0  (the  cylindrical  stagnation  point). 

CONDITIONS  FOR  ZERO  SECONDARY  FLOW 

With  secondary  flow  defined  as  flow  at  right  angles  to  the  invis- 
cid  streamlines,  we  see  that  it  is  given  by  the  function  G 7  of  Eqs.  (3) 
and  (4).  (Recall  that  these  are  written  for  "intrinsic”  coordinates, 
with  5  measured  along,  and  T)  normal  to,  the  inviscid  streamlines.) 

In  Eq.  (4),  every  term  contains  G  or  a  derivative  of  G,  except  for 
one  term,  proportioral  to  X,  and  another  proportional  to  u)^*  We  rule 
out  the  second  by  postulating  ze~o  spin.  In  order  to  make  X  vanish, 
we  first  recall  that 


Boundary  Layer  and  Flow  Control,  G.  V.  Lachmann,  ed.,  Pergamon 
Press,  1961  (hereafter  cited  in  the  text  as  B.L.F.C.). 
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d(£n  h^)  ^  dh^ 

|  x  5  T&n o  '  Vi^aT  ’  h2%i 

where  Is  the  geodesic  curvature  of  the  coordinate  line  of  constant 
TJ  (in  this  case  the  invi3cid  surface  streamline).  Geodesic  curvature 
of  a  surface  curve  at  a  point  P  is  the  curvature  of  the  projection  of 
that  curve  on  the  plane  tangent  to  the  surface  at  P.  Surface  curves 
are  called  geodesics  if  they  have  no  geodesic  curvature  and  hence  ap- 
*  pear  locally  straight  when  viewed  along  a  normal  to  the  surface  at  P. 

f  There  are  an  infinite  number  of  geodesic  curves  through  P,  each  with 

its  cwn  tangent  direction.  (A  familiar  example  is  given  by  the  "great 
circle"  routes  over  the  earth's  surface  from  one  city  to  various  others.) 

Hence,  if  =*  0  and  if  the  inviscid  flow  surface  streamlines  are 
geodesic  curves  of  the  body  (as  they  are,  for  example,  in  the  Newtonian 
hypersonic  flow  theory)  the  equation  for  G '  will  contain  no  term  inde¬ 
pendent  of  G Since  G ’  is  subject  to  homogeneous  boundary  conditions, 

G  r  «  0  is  the  indicated  solution  and  there  will  be  no  secondary  flow, 

SIMILAR  SOLUTIONS  FOR  SIMPLE  SPINNING  FLOWS 

When  (jj^  ^  0  it  is  still  possible  to  find  some  exact  similar  bound¬ 
ary-layer  solutions,  but  we  omit  discussion  of  these  here,  for  lack  of 

40 

time  and  space.  The  reader  is  referred  particularly  to  Moore  (1956) 
and  Mager  (1963)  for  excellent  reviews. 


VI.  CALCULATION  SCHEMES  FOR  ARBITRARY  30DY  SHAPES 

In  problems  in  which  the  right-hand  sides  of  £qs.  (3)  and  (4)  do 
not  vanish,  we  must  deal  with  partial  differential  equations  in  three 
independent  variables.  The  usual  approach  is  to  replace  these  equa¬ 
tions  by  a  sequence  of  ordinary  differential  equations  (with  indepen¬ 
dent  variable  z)  by  either  (1)  postulating  convergent  expansions  of  f' 
and  G '  as  power  series  in  x  and  y,  with  coefficient  functions  f/^(z), 
G^j(z)  determined  by  ordinary  differential  equations  (a  generalization 
of  tne  Blasius  series  approach  to  two-dimensional  problems),  or  (2)  ap¬ 
proximating  x-  and  y-derivatives  by  finite  differences,  and  solving 
ordinary  differential  equations  at  each  mesh  point  (x^,  y^),  (a  gen¬ 
eralization  of  Smith  and  Clutter's  approach  to  two-dimensional  prob¬ 
lems).  In  either  procedure,  careful  attention  must  be  paid  to  the  con¬ 
cepts  of  region  of  influence  and  domain  of  dependence  in  order  to  fa¬ 
cilitate  the  organization  of  computations  and  to  avoid  fundamental  er- 


SERIES  EXPANSIONS 

I  am  not  aware  of  any  published  study  of  the  application  of  the 
Blasius  series  method  to  general  three-dimensional  boundary- layer  prob¬ 
lems.  There  have  been  several  applications  to  the  computation  of  cross- 
flow  velocity  profiles  in  cases  which  are  governed  by  the  independence 
principle,  or  which  are  reduced  to  quasi-two-dimensional  problems  by 
assumption  of  weak  cross  flow.  Some  of  these  are  reviewed  by  Mager. 

It  appears  in  some  of  these  works  that  more  terms  in  the  x-expansions 
of  Blasius  are  needed  to  provide  satisfactorily  accurate  cross-flow 
profiles  than  are  needed  to  get  good  strearawise  profiles.  It  might  be 
worthwhile  to  look  into  the  possibilities  of  applying  this  method's 
self-similar  nodal  attachment  point  solutions  to  some  body  such  as  an 
ellipsoid. 


FIKITE-DIFFEREHSE  PROCEDURES 


Raetz  has  been  the  principal  contributor  to  the  technique  of  di¬ 
rect  numerical  assault  on  the  three-dimensional  laminar  boundary- layer 
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equations.  He  and  his  co-workers  have  computed  the  growth  and  stability 
of  laminar  boundary  layers  on  airplanes,  accounting  for  compressibility 
and  distributed  suction.  We  shall  not  describe  the  details  of  his  meth¬ 
od,  but  only  sketch  here  a  generalization  of  the  approach  applied  suc¬ 
cessfully  to  nonsimilar  two-dimensional  boundary  layers  by  Smith  and 
Clutter  (1962  and  1964).  This  method  employs  the  boundary- layer  equa¬ 
tions  in  the  same  transformation  as  we  use  here,  and  seems  a  little 
easier  to  explain.  Given  the  difficulties  of  guaranteeing  the  accura¬ 
cy  of  finite-difference  solutions  of  partial  differential  equations, 
it  might  actually  be  worthwhile  to  develop  this  procedure  into  a  work¬ 
ing  alternative  to  Raetz's  method,  and  to  compare  the  two  by  applying 
them  to  identical  problems. 

For  two-dimensional  or  axisymmetric  incompressible  flow.  Smith  and 
Clutter  (1963a)  integrated  the  equation 

f"+  (ai-i  +  h)  if  +  »(1  -  t -2)  -  *  (f '  tr  -  £'  i) 

for  specified  distributions  of  m  and  k  versus  x.  The  calculation 
started  with  initial  data  given  by  the  appropriate  self-similar  solu¬ 
tion  at  x  a  0,  and  thereafter  approximated  df  '/dx  and  df/dx  by  trail¬ 
ing  finite  differences.  At  each  discrete  value  of  x^,  an  ordinary  dif¬ 
ferential  equation  for  f(z,xn)  (in  which  the  values  of  f,  f'  etc.  at 
xn  ^  appear  as  variable  coefficients)  is  integrated  by  numerical  meth¬ 
ods  familiar  from  the  study  of  similar  boundary  layers.  This  consti¬ 
tutes  a  computationally  stable,  implicit  finite-difference  scheme,  the 
accuracy  of  which  could  be  controlled  by  restricting  the  permissible 

values  of  x  /Ax  .  In  their  second  paper,  the  method  is  extended  to 
n  n 

compressible  flow,  requiring  the  simultaneous  numerical  solution  of 
two  coupled  ordinary  differential  equations  at  each  x-value. 

The  method  suggested  here  is  analogous  to  the  former  in  that  fi¬ 
nite-difference  operations  would  be  used  for  both  d/dx  and  d/dy,  and 
a  coupled  pair  of  ordinary  differential  equations  integrated  at  each 
(x,y)  mesh  point.  Because  the  finite-differencing  must  not  be  done  in 
two  dimensions,  even  more  care  and  study  would  probably  be  required 
to  insure  computational  stability. 


t 
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The  finite-difference  procedure  would  be  started  at  an  "initial 
data  surface"  surrounding  a  nodal  attachment  point  or  a  sharp  edge  of 
attachment,  and  would  probably  employ  intrinsic  surface  coordinates 
[Eqs.  (3)  and  (4)1- 

One  point  which  might  require  special  care  in  the  finite-differ¬ 
encing  with  respect  to  y  is  related  to  the  fact  that  G '  can  change  sign 
in  the  boundary  layer.  Since  the  sign  of  G'  indicates  the  local  direc¬ 
tion  of  convection  by  the  secondary  flow,  and  since  a  fluid  particle 
carries  its  initial  data  downstream,  it  would  probably  be  necessary  to 
devise  a  differencing  scheme  which  uses  the  data  from  the  side  from 
which  the  flow  approaches  the  point  in  question.  If  the  inviscid  flow 
is  irrotational,  we  can  conveniently  let  the  Q  »  constant  curves  be 
equipctentials,  specifically  setting  §  =  §(£  =*  0)/U  ,  where  U  is  a 

CD  03 

constant  reference  speed.  Then  h  =  U  /U  and  \  +  n  =  0.  The  determi- 

1  03 

nation  of  h  is  then  directly  coupled  to  the  preliminary  task  of  com- 
^  * 

outing  the  inviscid  flow  over  the  surface.  The  coordinate  curves  of 
constant  (the  inviscid  streamlines)  might  rationally  be  spaced  (and 
h^  thereby  defined)  so  that 

i+1 

Uh2  dV, 
i 

is  independent  of  i,  when  the  integral  is  evaluated  along  an  equipoten- 
tlal.  This  subdivides  the  inviscid  surface  flow  which  "issues  from 
the  attachment  point"  into  equal  amounts  for  equal  increment  of  Tj. 

For  later  reference  we  note  that  in  this  coordinate  system,  h  measures 
the  lateral  spreading  of  inviscid  surface  streamlines --a  three-dimen¬ 
sional  effect  which  can  be  artificially  eliminated  by  a  generalized 
Mangier  transformation.  The  combination  m  +  h  measures  the  surface 
divergence  of  the  inviscid  velocity  field,  and  vanishes  for  problems 
in  which  the  inviscid  flow  is  two-dimensional  in  planes  parallel  to 
the  (necessarily  plane)  wall  on  which  the  boundary  layer  grows.  Such 
might  be  the  case  on  the  floor  of  a  wind-tunnel  turning  section. 

an  important  modern  reference  on  this  subject  is  Heaa  and  Smith 
(1963). 


f 


* 
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WEAK  CROSS-FLOW  METHODS:  PREVALENCE  PRINCIPLE 

A  great  many  approximate  calculations  of  three-dimensional  bound¬ 
ary  layers  take  advantage  of  the  fact  that  the  inviscid  streamlines 
may  not  deviate  too  much  from  surface  geodesics,  particularly  in  re¬ 
gions  of  favorable  streamwise  pressure  gradient,  and  of  course,  in  the 
vicinity  of  special  lines  of  body  symmetry.  In  such  cases  secondary 
flows  are  weak  and  have  small  effect  on  the  primary  or  streamwise  ve¬ 
locity  profiles.  This  is  called  the  "principe  de  prevalence"  by  Eichel- 
brenner  (1957),  and  the  corresponding  theory  might  be  more  prosaically 
termed  the  ‘‘weak  cross-flow  theory.'* 

In  its  "zeroth  approximation,"  the  weak  cross-flow  theory  sets  G 1 
and  G  equal  to  zero  in  £q.  (3),  which  thus  is  returned  almost  to  the 
form  treated  by  Smith  and  Clutter  (1963a),  namely 


f"  + 


f  m  +  1 


\ 

K) 


ff' 


hi 


f ,2) 


X_  (r;  I  df  7  _  f  P  5f  ^ 

hj  \  dx  / 


This  equation  can  be  integrated  along  each  inviscid  streamline.  The 
first  nonzero  estimate  of  secondary  velocities  is  then  obtained  from 
Ea.  (4),  in  which  only  terms  linear  in  G  and  its  derivatives  are  re¬ 
tained.  This  is  again  an  equation  which  can  be  integrated  along  in¬ 
viscid  streamlines  by  Smith  and  Clutter's  method. 


G"'+ 


+  *)«'  - 


f /2) 


+  2  <~J  x  (1  -  fO 


in  which  f  and  f '  are  taken  to  be  the  zeroth-order  streamwise  solu¬ 
tions.  The  y-derivativea  which  would  appear  in  the  equation  for  the 
next  approximation  to  f  are  then  approximated  by  finite  differences 
constructed  as  known  functions  of  x  by  use  of  the  lower  approximations 
to  f  and  G  on  neighboring  streamlines  (curves  of  constant  y),  and  once 
again  v»c  obtain  an  equation  which  can  be  integrated,  in  the  manner  of 
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Smith  and  Clutter,  alorg  the  invlscid  streamlines.  An  analogous  pro¬ 
cedure  can  be  applied  to  the  momentum-integral  equations,  to  convert 
them  from  partial  differential  equat: one  in  x  and  y  to  a  sequence  of 
approximate  ordinary  differential  equations  in  x.  This  procedure  is 
discussed  by  Eichelbrenner  and  by  Mager  (1964). 
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VII.  VORTICITY  OF  THREE-DIMENSIONAL  BOUNDARY  LAYERS 


Sometimes  we  can  establish  an  intuitive  grasp  of  complicated  flow 
phenomena  by  considering  the  sources  and  history  of  vorticity,  parti¬ 
cularly  in  the  case  c1  three-dimensional  flows.  Our  building  blocks 
are  the  vorticity  equation  and  certain  results  about  sources  of  vorti¬ 
city  at  solid  walls,  and  our  end  product-  is  an  understanding  of  the 
origins  and  development  of  secondary  vorticity  and  of  the  pattern  of 
vortex  lines  and  skin  friction  lines  on  the  wall. 

From  the  vorticity  equation  we  learn  that  there  is  only  one  ini¬ 
tial  volume  source  of  vorticity  in  an  initially  irrotational  flow  un¬ 
der  the  action  of  body  forces  derivable  from  a  potential.  This  source 
proportional  to  grad  p  X  grad  p  is  called  baroclinicity  by  meteorolo¬ 
gists,  and  it  vanishes  in  a  fluid  of  constant  density.  In  its  absence, 
vorticity  must  initially  appear  at  a  boundary  of  the  fluid  region. 

Once  in  the  fluid,  it  may  be  diffused  from  particle  to  particle  by 
viscous  torques,  the  direction  of  diffusion  being  along  the  normal  to 
the  wall  in  boundary- layer  flows.  Aside  from  the  effects  of  barocli¬ 
nicity  and  viscous  torques,  vorticity  is  simply  carried  along  with  the 
moving  fluid  in  such  a  way  as  to  conserve  the  flux  of  0  +  2<i  across 
any  surface  element  of  fixed  material  identity.  Many  implications  of 
this  are  exhibited  in  the  films  ’’Vorticity"  and  "Secondary  Flow"  in 
the  series  sponsored  by  the  National  Committee  for  Fluid  Mechanics 
Films. 

THE  SURFACE  SOURCE  OF  VORTICITY  IN  INCOMPRESSIBLE  FLOW 

The  diffusive  flux  of  vorticity  Q  across  a  surface  in  incompres- 
sible  flow  can  be  described  by  the  tensor  v  grad  fi,  in  the  sense  that 
the  vector  -n(v  grad  0)  gives  the  race  that  vorticity  is  diffused 
across  a  surface  normal  to  the  unit  vector  n,  into  the  region  into 
which  n  points.  If  the  surface  happens  to  be  a  solid  wall,  it  is 
natural  to  speak  of  -n(v  grad  Cl)  as  the  vorticity  source  strength  at 
a  point  on  the  wall.  Under  the  assumption  of  zero  slip,  and  with  the 
assumptions  of  boundary- layer  theory,  specifically  that  h..  =  h, C§,T1) , 

h2  “  3nd  h3  “  We  find 
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The  surface  source  of  tangential  components  of  vcrticity  is  pre¬ 
determined  by  the  inviscid  flow  over  the  body.  If  we  imagine  the  func¬ 
tion  P(§,T|)  to  define  a  pressure  or  potential  "hill,"  tangential  vor- 
ticity  is  fed  into  the  boundary  layer  at  a  rate  proportional  to  the 
slope  of  this  hill,  and  with  the  sense  of  rotation  which  would  be  taken 
by  a  ball  released  from  rest  to  start  rolling  down  the  hill. 

There  may  also  be  a  surface  source  of  normal  vorticity,  but  this 
is  not  predetermined  by  the  inviscid  flow.  Like  the  normal  component 
of  velocity  at  the  outer  edge  of  a  boundary  layer,  or  like  the  surface 
values  of  tangential  vorticity,  it  is  part  of  the  response  of  the 
boundary  layer  to  imposed  conditions- 

If  we  define  "primary"  or  "lateral"  vorticity  as  that  component 
of  0  which  is  tangential  to  the  wall  and  normal  to  the  surface  stream¬ 
lines  of  inviscid  flow,  and  "secondary"  or  "streamwise"  vorticity  as 
that  component  parallel  to  the  inviscid  surface  flow,  and  adopt  in¬ 
trinsic  surface  coordinates  in  which  V  =  0,  we  find  that  the  wall 

source  strength  of  primary  vorticity  is  -(U/hj)  dU/d§  and  that  the 

o 

wall  source  strength  of  secondary  vorticity  is  Hj,U  -  where 

is  the  geodesic  curvature  of  the  inviscid  surface  streamlines. 

Of  course,  these  expressions  give  only  the  tangential  components 
of  acceleration  of  the  inviscid  flow,  along  and  normal  to  its  surface 
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streamlines.  If  the  flow  were  unsteady,  there  would  be  a  wall  source 
of  primary  vorticity  proportional  to  -SU/St- 

Lighthill,  in  L.B.L. ,  pp.  83-84,  gives  an  interesting  heuristic 
derivation  of  these  results  without  making  reference  to  the  pressure 
field,  thereby  making  clear  that,  with  p  =  constant,  it  is  the  accel¬ 
eration  of  the  inviscid  flow  that  is  responsible  for  the  wall  sources 
of  vorticity. 

IKS  SURFACE  VORTICITY 

What  we  really  wish  to  know  in  a  typical  boundary- layer  problem 
is  the  distribution  of  surface  vorticity  which  results  from  the  pre¬ 
determined  distribution  of  wall  sources.  According  to  our  discussion 
in  Section  III,  we  can  see  that  the  surface  vorticity  at  a  point 
P(§»'n»0)  must  depend  upon  the  sources  distributed  throughout  its  "up¬ 
stream-opening"  region  of  influence.  This  region  is  defined  by  the 
normal  projection  on  the  wall  of  all  the  upstream  trajectories  of  par¬ 
ticles  which  arrive  in  the  boundary  layer  over  ?.  Actually,  not  all 
such  sources  have  equal  influence  on  the  surface  vorticity  at  P,  since 

the  contribution  of  a  given  upstream  source  is  proportional  to  its 

-1/2 

strength  and  to  a  damping  factor  (vt)  ,  where  t  Is  an  "average  time 
of  convection"  between  the  source  point  and  P. 

SURFACE  VORTKX  LINES  AND  SKIN  FRICTION  LIKES 

We  can  imagine  the  body  surface  to  be  covered  with  a  set  of  curves 
everywhere  tangent  to  the  surface  vorticity  vector.  Since  on  a  sta¬ 
tionary  surface  with  zero  slip,  the  Navier-Stokes  stress  relationship 
gives  a  surface  shear  stress  (skin  friction)  vector  related  to  the  sur¬ 
face  vorticity  by 


w 


n£ 


nB 


X  n 


we  see  that  the  family  of  skin  friction  lines,  which  are  everywhere 
tangent  to  t  ,  form  an  orthogonal  mesh  with  the  surface  vortex  lines. 
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The  skin  friction  lines  are  also  frequently  referred  to  as  "sur¬ 
face  streamlines,"  since  they  have,  according  to  the  boundary- layer 
approximations,  a  tangent  which  coincides  with  the  limiting  direction 
of  boundary- layer  streamlines,  as  Q  -  0.  That  Is 


lim 

C-0 


xiLat-Q 


0 

o 


but  the  indeterminacy  can  be  resolved,  except  at  isolated  singular 
points,  by  L'Hospital's  Rule,  to  give 


lim 

C-o 


v(I.Tl,C) 

ud^.O 


(av'*c)d,TUQ)  _  li 

du/aod.-H.o)  "  j1 


By  use  of  div  Q  =  0  and 


t  =  X 
w 


n),.  n  we  can  show  that 


so  that  the  wall  source  of  normal  vorticity  is  proportional  to  the 
circulation  of  the  skin  friction  field  around  the  point  of  observa¬ 
tion. 
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VIII.  TOPOLOGY  OF  THE  SKIN  FRICTION  LINES 


The  skin  friction  lines  of  a  steady  flow,  like  streamlines  in  the 
body  of  the  fluid,  can  run  together  or  cross  one  another  only  at  iso¬ 
lated  singular  points,  which  arise  where  the  vector  skin  friction  van¬ 
ishes  so  that  is  in  turn  indeterminate.  Such  "stagnation  points" 

of  the  skin  friction  field  are  points  of  flow  attachment  or  separation 
if  at  such  a  point  the  surface  divergence 

h,h2  (af  (lYV  +  2^  (hlV) 

is  positive  or  negative,  respectively.  Such  points  may  also  be  singu¬ 
lar  for  the  surface  vorticity  field,  in  which  "sources"  correspond  to 
points  of  attachment  of  vortex  filaments. 

The  topology  of  the  net  of  skin  friction  lines  and  vortex  lines 
is  dominated  by  these  singular  points,  and  special  studies  of  the  flow 
neighboring  such  a  point  have  been  made  by  R.  Legendre  (1955  and  1956) 
and  by  K.  Oswatir.sch  (1957).  We  give  here  a  very  brief  resume  of  the 
article  by  Oswatitsch,  which  contains  many  perspective  drawings  which 
greatly  facilitate  understanding. 

For  a  local  study  near  a  surface  point,  we  place  the  origin  of 
coordinates  at  that  point  (§  =  =  £  =  0)  and  adjust  the  (arbitrary) 

spacing  of  the  §-  and  H-coordinate  surfaces  so  that  h^(0,0,0)  = 
h?(0,0,0)  =  1.  Taylor  series  to  the  second  order  in  §,  H,  and  £  are 
postulated  for  u,  v,  and  w.  These  are  indicated  by  the  single  vector 
equation 


Many  English  authors  have  explored  the  possibilities  of  expan¬ 
sion  ether  than  the  simple  Taylor  series  in  the  vicinity  of  separa¬ 
tion.  For  a  recent  example  of  this  line  of  thought,  see  Brown  (1965). 


u(§,-\C>  =  u(0,0,u)  +  lf|  +  "|S+  Cf7 


+ 1  (c2  is.  +  J  is.  +  c2 

2  V  952  3"2  5C2' 

+  +  "C  +  C5  |^)  +  ... 

where  all  the  derivatives  of  u  are  evaluated  at  the  origin  (0,0,0). 
The  wall  is  stationary  and  impervious,  so  ud^.O)  =  0.  Thus,  at 
points  on  the  wall 
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From  the  continuity  equation,  we  next  find  that,  on  the  wall 

3w  _  d^w  _  d^w  _  n 

ac  asac  ' 

Finally,  we  differentiate  the  continuity  equation  with  respect  to  Q 
and  find  that,  on  the  wall 
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where  m. are  the  geodesic  curvatures  of  the  surface  coordinate 


The  momentum  equation  can  be  used  to  get 
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for  points  on  the  wall,  where  P  =  [p  +  p(A  -  1/2  a<  R  )^>  Finally,  we 
introduce  the  skin  friction  vector 


T 

W 


u  ||  (M,o) 


and  specify  that  the  origin  (0,0,0)  is  a  singular  point  at  which  =  0. 
All  this  leads  to  the  scalar  expressions 


nu(§»n,o  - 


uv(§,n,c)  = 


+ 

2  o§ 


TIC 


2  an 


or 

3ti 

w + 

aT  +  "  ‘ 

ar? 

_ 2  A 

an  ' 

«5T+  •*' 

nv(S,n,c) 


in  which  all  the  differential  coefficients  are  evaluated  at  (0,0,0). 

From  these  equations  follow  the  differential  equations  of  the 
streamlines  in  the  vicinity  of  (0,0,0),  namely. 
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Oswatitach  gives  Che  general  analytic  solutions  for  these  equations, 
from  which  we  sample  only  a  few  useful  results.  The  first  of  these 
has  to  do  with  the  skin  friction  lines  themselves,  whose  trajectories 
(at  Q  =  0)  are  governed  by  the  nonlinear  equation 


M 

d5 


We  first  see  whether  any  skin  friction  lines  enter  the  origin  by  seek¬ 
ing  the  singular  straight-line  solutions 
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These  exist  provided  that 
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This  has  two  distinct  real  values  if 


where 


St2  StA 

sr  ■  sr/ 


-  4 


A2  -  4J  >  0 


A 


dll  ’ 


dr^  dT2  dTj  dt2 

af  stT  "  W  df 


j 


-36- 


These  parameters,  together  with 


completely  determine  the  skin  friction  field  near  the  singular  point. 
There  are  few  general  physical  restrictions  on  the  values  of  the  diver¬ 
gence  cr  the  curl  (-11),  which  may  develop  in  the  skin  friction 

field,  so  we  may  expect  to  encounter  in  practice  singular  points  of 
great  variety. 

SYMMETRIC  SINGULAR  POINTS 

The  skin  friction  field  may  be  curl-free  (0  =  0)  everywhere  in 
very  special  cases  (e-g.,  axisymmetric  flow,  potential  inviscid  flow 
with  geodesic  surface  streamlines),  along  lines  of  symmetry  on  more 
generally  shaped  bodies  (e-g.,  an  ellipsoid  of  revolution  at  angle  of 
attack),  or  at  least  at  points  of  attachment  of  externally  irrotation- 
ai  flow.  Singular  points  encountered  under  these  conditions  exhibit 
a  symmetry,  in  that  the  singular  trajectories  through  them  are  mutu¬ 
ally  perpendicular,  so  that  the  §  and  T|  axes  can  be  rotated  to  coin¬ 
cide  with  them-  The  symmetry  then  yields 


and  the  general  trajectory  equation  becomes 


£1  = 
d§ 


-  ■  *  (£)/&) 


To  be  consistent  with  Oswatitsch's  terminology,  we  can  call  c  the  con¬ 
vergence  of  the  singular  point.  As  shown  In  Figs.  2  and  3,  the  solu¬ 
tion  curves 
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all  converge  at  the  origin  when  c  >  0  (a  node),  while  ail  approaching 
trajectories  except  two  (along  the  §  and  "h  axes)  diverge  away  from  the 
origin  when  c  <  0  (a  saddle  point)-  Convergence,  defined  in  this  sense, 
is  unfortunately  not  the  opposite  of  divergence,  A,  but  is  related  to 
it  by 


OT 

A  =  35-d  +  c) 

From  any  given  trajectory  plot,  one  can  quickly  evaluate  c  as  the 
slope  of  a  skin  friction  curve  at  the  point  where  it  crosses  the  line 
H  =  §•  It  is  easy  to  tell  the  sign  of  A  at  a  node  by  looking  at  the 
arrows  on  the  skin  friction  lines.  At  a  saddle  point  we  can  derive 
the  rule  of  thumb:  A  >  0  (attachment)  occurs  when  the  approaching 
skin  friction  lines  are  more  closely  packed  (near  their  axis  of  sym¬ 
metry)  than  the  departing  skin  friction  lines. 


ANGLE  OF  THE  ATIACKItC  OR  SEPARATING  STREAMLINE 

For  the  symmetric  singular  points,  0  =  0,  it  is  easy  to  find  the 
slope,  d£/d§,  of  the.  singular  streamline  which  leaves  the  wall  at  the 
singular  point.  We  suppose  that  11  =  0  is  a  plane  of  symmetry,  so  that 
dP/fiT)  =  0  (hence  dC/dlj  *  e») . 

The  general  equation  is 
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and  the  singular,  straight-line  solution,  is 
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where 
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The  first  factor  in  this  expression  is  the  familiar  angle  of  se¬ 
paration  of  a  two-dimensional  boundary  layer,  with  specified  values  of 
5P/d£  (>  0  for  separation)  and  (<  0,  but  bearing  no  unique  re¬ 
lation  to  We  see  that  at  sytrenetric  three-dimensional  singular 

points  it  is  modified  by  the  factor  1  +  c,  which  is  >  1  for  nodes,  and 
<  1  for  saddles.  Note  that  at  attachment  points  of  an  irrotational 
inviscid  flow,  BP/d§  =  0,  and  the  attaching  streamline  is  normal  to 
the  wall. 


MORE  GENERAL  SINGUIAR  POINTS 

As  illustrated  by  Legendre  and  Oswatitsch,  J  >  0  generally  denotes 

2 

a  node  and  J  <  0  a  saddle  point.  When  0  <  J  <  £  /4>  the  friction  lines 

2 

behave  qualitatively  as  around  a  symmetric  node,  while  for  J  >  £  /4 

2  2 

they  spiral  into  the  origin-  (Replace  £  by  Q  and  friction  lines  by 

vortex  lines,  and  the  last  sentence  is  again  true.)  If  J  is  greater 
2  2 

than  both  £  /4  and  0  /4,  both  skin  friction  lines  and  vortex  lines 
spiral.  Evidently  spiral  detachment  (J  >  £  /4  and  £  <  0)  occurs  quite 
commonly  in  practice,  notably  on  delta  wings-  Lighthill's  discussion 
and  sketches  on  pp.  74-82  of  L. B.L.  are  very  helpful. 

NUMBER  OF  MODES  AND  SADDLES  ON  A  CLOSED,  SIMPLY  CONNECTED  BODY 

li^hthill  cites  (p.  76,  L. B.L. )  the  following  topological  law, 
that  the  number  of  nodal  points  must  exceed  the  number  of  saddle  points 
by  2.  He  says  "one  may  argue  that  the  infinity  of  skin  friction  lines 
on  the  surface  must  begin  and  end  somewhere,  which  indicates  that  there 
is  at  least  one  nodal  point  of  attachment  and  one  nodal  point  of  separ¬ 
ation.  If  there  are  two  nodal  points  of  attachment,  the  skin  friction 
lines  from  each  must  somewhere  run  into  one  another,  and  so  have  to 
divide  at  a  saddle  point." 


-40- 


IX.  SEPARATION 


SEPARATION  UffiS 

From  Che  point  of  view  of  topology  of  skin  friction  curves,  un¬ 
separated  flow  exists  when  the  entire  body  surface  is  covered  by  fric¬ 
tion  lines  which  originate  at  upstream  nodal  points  of  attachment,  and 
disappear  into  an  equal  number  of  downstream  nodes  of  separation.  In 
practice,  this  situation  is  atypical;  experiments  ordinarily  show  the 
existence  of  other  nodal  "reattachmenf1  points,  at  which  the  incoming 
streamline  approaches  not  from  the  free  stream,  but  from  the  interior 
of  a  wake  bubble.  The  skin  friction  lines  from  the  upstream  attach- 
ment  node  and  those  from  the  wake  reattachment  node  run  towards  each 
other  and  hence  must  divide  at  a  saddle  point,  which  will  ordinarily 
be  a  saddle  point  of  separation.  The  singular  skin  friction  lines  ap¬ 
proaching  this  saddle  point  come  from  the  nodes  of  attachment  and  re¬ 
attachment.  The  singular  skin  friction  lines  departing  from  the  sad¬ 
dle  point  run  eventually  into  a  node  or  two  nodes  of  separation. 

The  latter  singular  friction  lines  separate  the  body  surface  in¬ 
to  a  region  which  is  wetted  by  the  upstream  flow,  and  a  region  which 
is  wetted  by  the  reattaching  wake  flow.  They  are  thus  logically  de¬ 
fined  as  the  separation  lines. 

The  concepts  of  a  separation  line  and  of  nodal  points  of  attach- 
+  ment  and  of  reattachraent,  a  considerable  variety  of  symmetric  singular 

points,  and  the  rule  regarding  the  relative  number  of  nodes  and  saddles 
are  all  illustrated  in  Fig.  4.  In  it  we  view  a  "foot-long  hot  dog," 
sagging  somewhat  in  the  middle.  A  uniform  flow  approaches  along  the 
x  axis,  and  the  hot  dog  is  supposed  to  possess  mirror  symmetry  in  the 
planes  y  =  0  and  z  =  0.  We  view  it  along  a  line  that  is  inclined 
slightly  upstream  and  spanwise  from  the  z  axis,  so  that  we  see  the 
forward  attachment  points  and  the  interesting  details  on  one  end.  The 
sketch  is  only  intended  to  appear  plausible,  and  does  not  come  from 
an  actual  experiment  or  calculation.  In  particular,  if  the  hot  dog 
were  drooped  more  in  the  middle  and  less  at  the  ends,  many  more  sing¬ 
ular  points  might  be  expected  to  appear  in  those  regions.  (These 
guesses  derive  from  the  paper  on  saddle  points  of  attachment  by  Davey 
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(referred  to  in  Section  VI),  and  from  drawings  derived  from  observed 
flow  on  a  yawed  ellipsoid  by  Eichelbrenner ,  and  exhibited  in  the  sur¬ 
vey  article  by  Cooke  and  Brebner  in  B.L.F.C. 

INFLUENCE  Of  SEPARATION  OH  BOUNDARY -LAYER  COMP UTAH QIC 

Separation  limits  the  applicability  of  boundary- layer  theory  in 
two  ways.  The  first  can  be  anticipated  within  the  framework  of  the 
theory  and  our  previous  discussion  of  domains  of  dependence  and  appro¬ 
priately  chosen  initial  data  arcs.  From  these  considerations  it  is 
clear  that  we  cannot  "march'1  our  computations  across  a  separation  line, 
since  the  flow  on  the  downstream  side  of  such  a  line  falls  into  the  re¬ 
gion  cf  influence  of  the  unknown  reattaching  flow. 

The  next  question  is  whether  we  can  march  up  to  the  separation 
line  from  the  wake-wetted  side  as  well  as  from  the  upstream-wetted 
side,  and  thereby  compute  the  boundary-layer  development  over  the  en¬ 
tire  body.  The  answer  appears  to  be  almost  always  negative-  One  dif¬ 
ficulty  is  that  we  do  not  know,  a  priori,  where  the  wake-reattachment 
nodes  will  be.  They  do  not  in  general  coincide  with  any  singular  fea¬ 
tures  of  an  ii rotational  inviscid  solution  (and  in  practice  they  tend 
to  wander  temporally,  at  any  Reynolds  numbers  for  which  the  boundary 
layers  may  be  expected  to  be  thin).  Furthermore,  the  fluid  attaching 
at  these  points,  unlike  the  free  stream  fluid  attaching  at  upstream 
nodal  points,  has  had  a  recent  experience  of  viscous  action,  so  that 
we  cannot  ordinarily  expect  to  find  self-contained  similar  solutions 
for  the  vicinity  of  the  reattachraent  node.  Finally,  even  if  these 
difficulties  and  those  of  determining  the  inviscid  surface  velocities 
in  the  wake-wetted  region  could  be  surmounted,  there  remains  the  pos¬ 
sibility  that  near  the  separation  line  the  direction  of  dominant  vis¬ 
cous  diffusion  differs  significantly  from  that  of  the  wall  normal, 
thus  invalidating  our  boundary- layer  equations  locally. 

The  second  adverse  effect  of  separation  is  its  displacement  ef¬ 
fect,  which  causes  the  actual  pressure  distribution  to  deviate  from 
that  of  inviscid  theory  over  a  region  which  may  extend  significantly 
upstream  of  the  computed  line  cf  separation.  This  effect  would  then 
render  inaccurate  the  boundary- layer  computations  in  this  vicinity, 
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and  in  particular  their  prediction  of  the  location  of  the  separation 
line. 

Though  little  direct  evidence  as  yet  exists,  there  is  reason  to 
hope  that  this  latter  effect  may  be  less  damaging  in  many  types  of 
three-dimensional  flow  than  in  two-dimensional  flow.  These  intuitive 
hopes  derive  from  the  expectations  chat  (1)  the  angle  of  separation 
may  be  quite  small  at  and  near  the  saddle  points  of  separation  (A  be¬ 
ing  small  there) ;  (2)  the  line  of  separation  may  be  considerably  "swept 
back"  relative  to  the  local  inviscid  flow;  (3)  inviscid  flow  can  to 
some  extent  go  around  the  regions  of  vigorous  outflow  near  the  nodal 
separation  points. 

To  set  against  this  optimistic  view  is  the  observation  that  spi¬ 
ral  nodes  of  separation  are  often  situated  at  the  feet  of  tightly  wound 
and  vigorous  "trailing  vortices"  which  may  induce  profound  perturba¬ 
tions  of  the  neighboring  inviscid  fields. 
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X.  MOMENTUM-INTEGRAL  METHODS 


Granting  the  fact  that  the  direct  numerical  integration  of  Eos- 
(1)  and  (2)  or  of  Eqs.  (3)  and  (4)  will  be  very  lengthy  even  if  the 
method  is  computationally  stable  and  well-posed,  we  may  wish  to  inves¬ 
tigate  approximate  methods  which  satisfy  the  equations  in  a  (possibly 
weighted)  mean  sense  at  each  x  and  y.  These  are  the  momentum-integral 
methods,  of  which  a  particularly  comprehensive  review  is  given  by 
Lindfield,  Pinsent,  and  Pinsent  in  B.L.F.C.  We  mry  derive  these  equa¬ 
tions  directly  from  Eqs-  (1)  and  (2)  or  Eqs.  (3)  and  (4)  as  follows. 

THE  MOMENTUM-INTEGRAL  EQUATIONS  IN  INTRINSIC  COORDINATES 

Holding  x  and  v  constant,  we  integrate  Eqs.  (3)  and  (4)  over  z 
from  0  to  and  employ  the  appropriate  boundary  conditions  leer  f  and 
G-  For  an  impervious  wall,  we  have 

f(0)  =  f'(0)  =  G(0)  =  G'(0)  =  0 

The  matching  conditions  are  f  '(<»)  =  1  (hence  fff(°°)  =  0,  G/(»)  =  G/,(®) 

=  0).  Integration  is  assisted  by  the  device  of  setting 

f"  =  -(1  -  f ')  ' 


so  that  in 
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the  first  term  on  the  right  vanishes  at  both  limits. 

We  introduce  the  -  'Us 
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9ll(x’y)  s  J  (l-fOf'dz,  ei2(x,y)  s  J  (l-f')G'dz 
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For  easy  comparison  with  the  equation  of  Mager  (1364),  we  collect  terms 
in  0^  and  0^>  and  note  that 
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where  0^  is  the  (^-component  0f  the  inviscid-f low  vorticity.  This  al¬ 
lows  a  neat  combination  of  and  into  one  terra  which  vanishes  in 
many  applications.  Then  we  have 


i-  ?!u  +  J_  fu  .  m  1  (I3s.ii  +  2A  ..  +  e  ;  +  m,\ 

hj^  dx  n2  dy  x  h  x  tv  2  /II  11  22  1/ 


b  {(Is  +  2X>12}  +  <n3  +  2“3)  I  62  (5) 
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and  by  a  similar  calculation  on  the  cross-stream  momentum  equation 


1_ 
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-  {(?  +  2X)  °22  -  K8u 


-  <"3  +  2”3>  3  61 
+  e22)  +  nsj  (6) 


Here  we  have  two  coupled  firat-order  partial  differential  equa¬ 
tions  for  the  six  integral  quantities  (6^,  62,  8^,  ®21*  anc*  Q22) 

and  the  two  wall-derivatives  f"(0)  and  G*(G). 


COMPATIBILITY  CONDITIONS 

We  have  also  at  our  disposal  a  variety  of  algebraic  "compatibil¬ 
ity  conditions"  obtained  from  Eqs.  (3)  and  (4)  and  from  the  z-deriva- 
tives  of  these  equations,  evaluated  at  the  wall  (z  =  0) : 

First  Compatibility  Condition 


f'"(0) 


-  and  G"'(0) 

hl 


ID- Lx 

2  - 

u 


Second  Compatibility  Condition 


f'tO)  -  0  and  G""  -  xf"(0) 

(In  Lindfield,  Pinsent,  and  Pinsent  these  conditions  and  the  momentum- 
integral  equations  are  generalized  to  include  the  effects  of  wall  suc¬ 
tion  or  blowing.) 

SINGLY  INFINITE  FAMILIES  OF  VELOCITY  PROFILES 


We  follow  Head  (1961)  (whose  excellent  survey  of  integral  methods 
for  two-dimensional  flows  precedes  the  article  by  Lindfield  et  al.  in 
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Vol.  2  of  B.L.F.C.)  in  calling  a  family  of  velocity  profiles  singly 
infinite  if  the  shape  of  a  member  profile,  with  respect  to  a  linearly 
stretched  z  coordinate,  is  determined  by  a  single  parameter  or  shape 
factor.  The  profile,  expressed  in  terms  of  z  itself,  then  requires 
two  parameters,  the  shape  factor  and  the  scale  factor. 

Mathematically,  we  write 


A(x,y) 


Before  beginning  a  detailed  discussion  of  solution  procedures  and  ap 
proximating  expressions  foi  f/  and  G  ,,  we  note  the  following  gain  in 
simplicity  from  £qs.  (3)  and  (4)  to  Eqs-  (5)  and  (6). 

Our  proposed  solution  method  for  Eqs-  (3)  and  (4)  was  to  employ 
finite  differences  in  x  and  y,  to  obtain  two  coupled  ordinary  differ¬ 
ential  equations  at  each  (x,y)  mesh  point-  These  latter  were  to  be 
solved  numerically,  subject  to  the  split  boundary  conditions  (at  z  =  0 
and  z  =  - 

If  v/e  repeat  this  procedure  for  the  set  (5)  and  (6),  we  have,  af¬ 
ter  finite-differencing  in  x  and  y,  two  coupled  algebraic  equations  to 
be  solved  at  each  mesh  point-  Clearly  the  computational  savings  are 
great,  even  though  the  computational  effort  of  the  momentum-integral 
method  is  considerably  greater  here  than  in  two-dimensional  problems- 


WEIGHTED  INTEGRAL  EQUATION 

For  an  alternative  to  the  compatibility  conditions  we  may  use  an¬ 
other  pair  of  differential  equations  to  relate  the  shape  factors  to  the 
scale  factors.  Such  equations  may  be  constructed  by  first  multiplying 
the  momentum  equations  (3)  and  (4)  by  a  '‘weighting"  factor  (or  factors) 
and  then  integrating  over  z-  These  new  partial  differential  equations 
will  be  coupled  to  the  unweighted  momentum- integral  equations,  and  the 
substitution  of  them  for  the  compatibility  conditions  will  Improve  the 
accuracy  of  the  momentum- integral  method  only  in  return  for  an  increased 
computational  effort.  Head  gives  a  number  of  examples,  comparing  the 
accuracy  of  various  alternative  methods  for  some  important  and  particu¬ 
larly  trying  two-dimensional  boundary- layer  problems- 

One's  hope  for  a  gain  in  accuracy  through  introduction  of  the 
weighted  momentum-integral  equations  is  strongest  in  cases  involving 
a  sudden  change  in  pressure  gradient  (either  streamwise  or  cross¬ 
stream).  In  such  cases  the  first  compatibility  conditions,  which  al¬ 
ways  insist  on  an  exactly  correct  voiticity  source  strength  at  the 
wall,  may  wrench  the  rest  of  the  profile  severely  out  of  shape  in  the 
region  just  downstream  of  the  sudden  change-  This  effect  is  worst 
when  we  employ  singly  Infinite  profile  families,  which  are  incapable 
of  the  sort  of  shape  adjustments  which  occur  in  the  real  boundary  layer 


-49- 


t 

2 


1 


as  the  new  vorticity  gradually  diffuses  cue  from  Che  wall.  The  typi¬ 
cal  symptom  of  inaccuracy  is  a  discontinuity  or  excessively  rapid  change 
ci  surface  vorticity  (skin  friction)  and  displacement  thickness  at  the 
location  of  the  sudden  change  of  pressure  gradient.  In  these  regions 
the  curve  of  vorticity  versus  z  has  the  correct  slope,  but  a  very  in¬ 
accurate  intercept,  at  z  =  0- 

The  weighted  momentum- integral  equations  ordinarily  give  special 
emphasis  to  the  role  of  convection  in  the  boundary  layer,  the  usual 
weighting  factors  being  f'  for  Eq.  (3)  and  G'  for  Eq.  (4).  In  a  re¬ 
gion  of  sudden  pressure  change  these  equations  require  a  good  descrip¬ 
tion  of  the  quasi-inviscid  response  of  velocities  to  pressure  gradi¬ 
ents,  and  since  that  response  is  dominant  in  the  outer  regions  into 
which  the  fresh  vorticity  has  not  yet  diffused,  they  give  a  better  pic¬ 
ture  of  the  response  of  "overall"  boundary- layer  profile  shape  and 
thickness.  Close  to  the  wall  the  vorticity  is  forced  to  assume  reason¬ 
ably  correct  values  since  in  a  singly  infinite  profile  family  the  wall 
vorticity  and  the  boundary  layer  thickness  are  usually  intimately  con¬ 
nected.  The  slope  of  the  vorticity  profile  at  z  =  0,  and  hence  the  vor¬ 
ticity  source  strength,  will  probably  be  given  poorly  by  this  method  in 
the  regions  of  sudden  adjustment,  but  unless  great  accuracy  is  needed 
in  the  velocity  profile  curvature,  this  error  may  be  quite  tolerable- 

SAMPLE  SIMPLY  INFINITE  PROFILE  FAMILIES 

Head  and  Lindfield  et  al.  present  good  discussions  of  the  way  in 
which  profile  families  can  be  generated  by  using  polynomials,  trans¬ 
cendental  functions,  or  numerical  functions  of  the  stretched  normal 
distance.  The  usual  procedure  is  to  make  a  linear  combination  of  two 
functions,  one  of  which  is  selected  to  give  a  fair  approximation  to 
the  profile  at  attachment  points  and  in  ti>2  neighboring  regions  in 
which  the  streamwise  pressure  gradient  is  favorable  and  in  which  the 
inviscid  streamline  curvature  has  the  same  sign  as  at  attachment.  Sam¬ 
ples  of  such  "first"  functions,  from  the  work  of  Timman  (1951)  and  Zaat 
(1956)  are 
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erf  <t>  - 


0  =- 


£ 
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which  satisfies 


f'CO)  =  0,  f'(0)  =  f'(0)  =  0 

f'(o)  =  1,  £'  =  f"'  =  f^n)(»)  =  0 

and 

-Y2 

G'(Y)  =  Ye 

*  *  s 

which  gives 

G/(0)  -  0,  G'(0)  =  j,  G''(0)  -  0 

G'(°°)  =  G"^  -  G^n)(®)  »  0 


These  functions,  which  fulfill  the  boundary  conditions  on  £'  and  G #, 
are  given  in  Fig.  5-  The  fact  that  f£' (0)  =  G^' (0)  =  0  has  an  impor¬ 
tant  bearing  on  the  algebraic  simplicity  of  the  subsequent  theory  when 
the  first  compatibility  equations  are  used. 

The  second  component  of  the  profile  must  now  vanish  at  z  =  0  and 
z  =  but  can  otherwise  be  chosen  arbitrarily  with  the  special  objec¬ 
tive  of  imitating  the  profound  changes  in  profile  shape  which  accompany 
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Fig.  5  —  Zaat  and  Tinman’s  profile  functions. 

the  onset  of  adverse  streanwise  pressure  gradients  and  changes  in  the 
sign  of  the  inviscid  streamline  curvature.  To  achieve  this  objective 
in  both  f '  and  G  \  Zaat  and  Tinman  propose  the  function 

*2  "  \  <er£c  0  -  e‘4  }  +  ^  0e"<S 
and  let  be  this  function  of  This  gives 

f2'<0)  -  0,  f-(0)  -  -  — p-,  f,,(0)  =  1 

f2'(»)  =  0,  f2  »  f2'  -  f2n)(»>  »  0 


Thus  we  could  write 
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f'(x>y>2)  =  £'  (|)  -  X*2  (f) 

and 

C'(x,y,2)  -  AoG'(f)-I  C'  (|) 

It  can  be  readily  verified  that  both  of  these  combinations  can  describe 
"two-sided''  as  well  as  "one-sided"  profiles,  as  a  result  of  variation 
in  x  and  £• 

The  extra  coefficient  Aq  is  permissible  as  far  as  satisfaction 
of  the  boundary  conditions  of  G'  is  concerned,  and  may  (I  think)  be 
employed  to  get  the  calculation  started  in  the  vicinity  of  an  attach¬ 
ment  point.  In  such  a  region  exact  cross-flow  profiles  are  given  in 
terms  of  the  similar  solutions  of  Howarth,  being 

G '  s 

where  x,  y,  c,  g 1  and  f 1  are  as  defined  in  Section  VI.  While  the  func¬ 
tions  g  1  and  f '  are  universal  in  the  neighborhood  of  the  attachment 
point,  we  see  that  G/  will  vary  with  x/y  along  the  initial  data  arc. 
Since  Zaat  and  Tinman's  G/  has  about  the  same  shape  as  gr  -  f we 

might  take  E  =  0  on  the  initial  data  arc,  and  let  A  **  (x  /cy  +  cy  /x  ) 
°  ’  o  o  Jo  Jo  o 

DOUBLY  INFINITE  PROFILE  FAMILIES 

While  the  singly  infinite  families  of  profiles  which  are  composed 
from  Zaat  and  Tinman's  functions  may  assume  at  least  qualitatively  rea¬ 
sonable  shapes  in  interesting  problems,  particularly  if  we  employ  the 
"gentler"  but  more  tedious  procedure  with  four  integral  equations  and 
ignore  the  compatibility  conditions,  they  do  not  have  enough  flexibil¬ 
ity  to  permit  accurate  t -presentation  of  velocity  profiles,  such  as  are 
needed  in  attempts  to  predict  separation  or  laminar  stability.  Head 
and  Lindfield  et  al.  show  how  this  situation  may  be  improved  by  the 


introduction  of  doubly  infinite  profile  families,  which  are  character¬ 
ized  by  two  shape  factors  and  a  scale  factor  for  each  family. 

If  two  new  shape  factors  are  introduced,  then  we  need  six  equa¬ 
tions  in  all.  These  may  be  the  two  momentum- integral  equacions,  two 
weighted  momentum-integral  equations  and  the  first  compatibility  con¬ 
ditions.  Obviously  other  combinations  can  be  used,  and  Head  has  tested 
some  of  these  for  two-dimensional  flows-  If  one  has  decided  to  use  the 
four  integral  equations  for  singly  infinite  profile  families,  not  too 
much  extra  work  is  involved  in  going  to  doubly  infinite  profile  fami¬ 
lies. 

In  fact,  for  three-dimensional  boundary  layers,  one  might  expect 
that  the  t'  and  G  '*  profiles  could  porsibly  “share"  some  parameters, 
so  that  doubly  infinite  families  could  be  used  without  the  need  for 
more  equations-  That  is,  we  might  be  able  to  extend  our  profiles  to 


G '  «  AG'  -  EG.'  -  XG,' 

O  1  4  J 


but  it  would  require  considerable  experience  to  discover  what,  if  any, 
expressions  for  f^  and  G^  might  make  this  scheme  profitable- 


SOLUTIONOF  THE  MOMENTUM-INTEGRAL  AND  ENERGY-INTEGRAL  EQUATIONS 


Even  though  the  momentum-integral  equations  can  presumably  be 
solved  numerically  with  much  less  effort  than  the  full  boundary- layer 
equations,  they  have  seldom  been  applied  to  truly  three-dimensional 
problems.  The  exceptional  examples  seem  to  be  discussed  by  Zaat 
(1956)  and  Eichelbrenner  (1957),  and  both  of  these  employed  solution 
methods  adapted  only  to  regions  with  weak  cross  flows,  Lindfield, 
Pinsent,  and  Pinsent  (1961)  give  a  general  discussion  of  the  use  of 
weighted  integral  equations,  but  work  no  examples  in  that  publics- 


tion.  Only  Tinman  (1951)  and  Hager  (1957)  appear  to  have  studied  the 
characteristics  of  the  partial  differential  Eqs.  (7)  and  (8),  and  no 
study  seems  to  have  been  made  of  the  characteristics  of  the  fourth- 
order  system  of  two  momentum- integral  plus  two  energy-integral  equa¬ 
tions-  I  have  made  various  attempts  to  fill  this  gap  in  this  Memo¬ 
randum,  but  conclude  that  the  topic  is  too  lengthy  and  coo  specific 
to  be  worth  including. 

Derivation  of  Weighted  Momentum-Integral  Equations 

If  we  weight  Eq.  (3)  by  f *  and  then  integrate  over  z,  the  terms 
on  the  right-hand  side  contribute  the  derivative  terms  in  the  stream- 
wise-  energy-  integral  equation.  These  are 

„  2  /!_  8en  .  l_  5si2\ 

2  yh^  dx  h2  by  J 

where 

en  3  f  f,<1  -  £,2> dz 

o 

CO 

e12  "  J  G'(i  '  f,Z)  dz 

o 


If  we  try,  following  Lindfield  et  al.,  to  weight  Eq.  (4)  also  by 
f we  find  it  is  Impossible  to  extract  the  differentiations  in  x  and 
y  from  under  the  integral  sign.  (I  can  only  conclude  that  Lindfield 
et  al.  succeeded  in  doing  this  by  virtue  of  an  error  in  calculus.) 

The  resulting  equation  would  thus  be  unbearably  tedious  to  employ. 
However,  if  we  weight  Eq.  (4)  by  the  velocity  G',  no  such  difficulty 
arises.  The  derivative  terms  are 
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where 


e22  "  -fG'3d2 

o 

The  momentum  and  energy  integral  equations  may  thus  be  summarized 
as 
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By  using  the  assumed  profile  families,  either  singly  or  doubly 
infinite,  we  can  express  all  quantities  in  these  equations  as  func¬ 
tions  of  the  two  scale  factors,  6  and  A.  and  the  shape  factors  (two 
for  singly  infinite  families;  four  for  doubly  infinite  families).  If 
we  employ  doubly  infinite  families,  we  assume  herewith  that  appropri¬ 
ate  compatibility  conditions  have  been  employed  to  eliminate  two  of 
the  shape  factors  in  favor  of  the  other  two  shape  factors  and  the  two 
scale  factors.  We  call  the  remaining  shape  factors  x  and  £.  Then, 
in  the  case  of  greatest  generality  (with  shared  shape  factors),  we 
would  know 
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=  81]L(S^,X,I),  9??  =  81?(f,A,x,D,  etc. 
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we  would  write 


11 


3x 


S911  36  .  5911  oA  S611  3^  5911  3E 

36  3x  T  3A  ox  3x  3x  3E  ox 


and  similar  expressions  for  the  ocher  derivatives.  The  resulting  four 
equations  for  6,  x>  and  2  can  he  abbreviated  by  use  of  an  index  no¬ 
tation  (summation  over  repeated  indices)  as 


3u.  3u 


a.  .  v — “•  +  b.  .  v — L  +  C.  -  0, 
ij  3x  xj  3y  i 


ft  =  1,  2,  3,  V 
Ij  =  1,  2,  3,  4J 


(9) 


We  identify 


ui  =  6- 


=5 


=  X 


%  -  E 


IJ 


1_ 

h. 


aeu 

30n 

fiii 

36u 

as 

3A 

dE 

3021 

3621 

S021 

as 

3A 

dX 

3Z 

Ssu 

3en, 

li. 

dell 

dell 

as 

d£ 

3X 

3E 

5e21 

3S21 

8e21 

0e21 

36 

3£ 

3X 

3E 

and 


-57- 


b.  .  =  — 


Difficulties  of  Further  Discussion 

One  can  easily  carry  on  formally,  writing  expressions  for  the 
characteristic  directions  of  the  set  Eq«  (9)  to  show  that  there  are 
either  0,  2,  or  4  real  characteristics  at  each  point,  and  to  indi¬ 
cate  plausibly  that  in  cases  of  weak  cross  flow  all  real  characteris¬ 
tic  slopes  are  close  to  that  of  the  inviscid  streamline.  The  essen¬ 
tial  arbitrariness  of  the  momentum-integral  methods  precludes  more 
precise  statements  for  the  c^se  of  general  cross-flow  magnitude.  This 

makes  the  numerical  values  of  the  coefficients  a. .  and  b. ,  a  result 

J-j 

not  only  of  physical  necessity  but  also  of  personal  choice  (exercised 
in  the  paramctrization  of  velocity  profiles).  Furthermore,  if  condi¬ 
tions  exist  in  which  some  of  the  characteristics  are  imaginary,  the 
criterion  for  such  conditions  does  not  seem  to  be  expressible  in  terms 
of  a  simple,  physical  concept  (such  as  that  of  subsonic  versus  super¬ 
sonic  flow  in  gas  dynamics).  A  feeling  for  the  behavior  of  the  char¬ 
acteristic  curves  appears  to  come  only  from  extensive  experience  with 
actual  calculations,  and  this  experience  does  not  yet  exist. 

Difficulties  of  Application  of  the  Method  of  Characteristics 

While  in  principle  integration  along  characteristics  (when  these 
are  real)  is  the  preferred  method  of  solution  of  equation  sets  such  as 
Eq.  (9),  numerical  execution  is  awkward  (leading  to  either  poor  accu¬ 
racy  or  an  excessive  amount  of  computation)  when  the  characteristic 
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lines  all  have  nearly  the  same  slope,  as  in  the  weak  cross-flow  case- 
Such  cases  bep  for  a  different  scheme  of  approximation-  Hager  (1957) 
amplifies  these  points  somewhat. 

Furthermore,  although  they  are  described  in  principle  by  Courant 
and  Friedrichs  (1948),  numerical  procedures  become  very  tedious  whe: 
as  many  as  four  real  characteristics  exist,  as  may  be  the  case  for  set 
(9)- 


XL.  EFFECTS  OP  COMPRESSIBILITY 


In  keeping  with  our  announced  intention  to  concern  ourselves  pri¬ 
marily  with  geometrical  rather  than  thermodynamical  complications  of 
boundary- layer  theory,  the  following  consents  will  be  very  brief.  A 
more  comprehensive  survey  of  the  formal  aspects  can  be  found  in  Hager's 
review  article. 

Compressibility  can  give  rise  to,  or  modify,  secondary  flows  in 
three  essential  ways.  One  is  seen  in  the  onset  of  differential  "buoy¬ 
ancy"  accelerations,  which  arise  in  a  variable-density  boundary  layer 
under  the  action  of  velocity-independent  force  fields  (such  as  centrif¬ 
ugal  or  gravity  fields).  A  sample  problem  would  be  to  determine  the 
boundary  layer  on  a  strongly  heated  body  that  is  moving  slowly  in  a 
horizontal  direction.  A  survey  of  laminar  flows  with  body  forces  is 
given  by  Ostracb  (1964)  in  the  same  volume  which  contains  Mager's  sur¬ 
vey.  A  second  is  easily  visualized  in  terms  of  the  simple  quasi-in- 

viscid  force-balance  explanation  of  streamline  curvature  in  a  three- 

2 

dimensional  boundary  layer.  Since  it  is  pu  times  the  geodesic  curva¬ 
ture  -  v.ich  (nearly)  balances  the  transverse  pressure,  gradient,  the  cur¬ 
vature  of  boundary- layer  streamlines  will  be  increased  by  processes 
such  as  viscous  dissipation,  which  reduce  p,  and  decreased  by  process¬ 
es  such  as  extreme  wall  cooling,  which  increase  p.  In  an  article  on 
the  applicability  of  Stewartson's  transformation  to  three-dimensional 
boundary  layers,  J.  C.  Cooke  (1961)  deduced  the  rule  of  thumb  that  if 
(a)  Pr  »  1,  (b)  n  «  T,  (c)  cross  flows  are  almost  negligible,  (d)  the 
wall  is  insulated,  and  (e)  the  external  flow  Mach  number  is  moderate, 
then  the  compressible-boundary- layer  flow  corresponds  to  an  incompres¬ 
sible  flow  over  the  same  surface,  but  with  all  pressure  gradients  en¬ 
hanced  by  the  factor  1  +  (y  -  1)M^ /2.  On  the  other  hand,  for  the 
boundary  layers  on  highly  cooled  wall*  of  yawed  blunted  cones  in  hy¬ 
personic  flow,  Vaglio-Laurin  (1959)  has  shown  that  a  weak-crots-f low 
assumption,  which  decouples  the  streawise-momentum  equation  from  the 
cross-flow  momentum  equation  and  allcws  the  former  to  be  reduced  to 

two-dimensional  form  by  Mangier' a  transformation,  is  often  justified 

2 

oy  the  near-constancy  of  pu  across  the  boundary  layer.  The  third  and 


least-studied  way  in  which  compressibility  may  engender  secondary  flows 
involves  the  displacement  effect  of  the  boundary  layer.  While  it  has 
been  possible  to  incorporate  this  effect  fairly  successfully  into  the 
theory  of  two-dimensional  hypersonic  boundary  layers,  the  only  litera¬ 
ture  dealing  with  three-dimensional  cases  (other  than  axisymmetric 
cases)  appear  to  refer  to  corner-flow,  or  boundary-region,  problems. 

While  we  cannot  elaborate  here  on  the  subject  of  displacement  ef¬ 
fects,  we  shall  state  the  definition  of  displacement  thickness  of  a 
three-dimensional  steady  boundary  layer,  following  Lighthill  (1958). 

We  refer  to  Fig.  6,  which  attempts  to  show  two  surface  streamlines  of 
the  inviscid  flow,  denoted  by  T|  and  T|  +  dl]  in  intrinsic  coordinates, 
and  originating  at  a  nodal  point  of  attachment,  0.  Ove.7  each  of  these 
streamlines  we  sketch  a  "wall  of  surface  normals,"  the  two  walls  meet¬ 
ing  over  0.  At  the  point  (|,1J)  we  erect  a  rectangle  with  base  h^  dll 
along  the  §  =  constant  curve,  and  height  6*  along  the  local  normal 
(C-axis). 


0 


Fig.  6  --  Control  volume  which  defines  displacement  thickness. 
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We  may  define  the  displacement  thickness,  6^,  at  §,  “  by  equating 
the  inviscid  mass  flow  through  the  rectangle,  p^Ul^  dT)  fcj.,  to  the  mass 
flow  "displaced  by  boundary- layer  action"  between  the  attachment  node 
and  §,  T).  Part  of  this  displaced  flow  is  accounted  for  by  the  deficit 
of  streamwise  flux  in  the  boundary  layer  at  §,  f).  This  deficit,  reck- 
oned  as  a  positive  number,  is  h^  dtj  jq  (p^U  -  Pu)  d£.  The  remainder 
of  the  displaced  flow  has  leaked  out  through  "sidewalls*'  of  normals, 
between  0  and  §,  due  to  boundary- layer  cross  flow.  Between  and 
+  d§'  this  deficit,  again  reckoned  as  a  positive  quantity,  is  d"H 
3/dT]  [h^  d£'  pv  d£j,  and  the  accumulated  deficit  between  boundary- 
layer  attachment  and  the  point  of  interest  is 


dTl  f  [h1(§,,Tl)  J  p(§')v(§')  dC]  dV 


Thus  we  get 


peU6*h2  d*  =  h2  d11  J  (PeU  “  Pu)  dC  +  J  fjf  j  Pv  dc]  &¥ 

lo  o  Jo“‘ 

Introducing  the  "displacement  integrals  of  the  primary  and  secon¬ 
dary  flows," 


we  get 


6i  B  J  (l  -  fs)  ^ 


62  »  J  f5"t 

A  *  /-> 
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6*  -  *i  *  «* 

In  closing  these  brief  remarks,  we  enter  the  observation  that 
transformation  methods,  intended  to  carry  compressible-flow  problems 
over  into  previously  solved  incompressible-flow  problems,  appear  to 
play  a  relatively  minor  role  in  the  theory  of  three-dimensional  bound¬ 
ary  layers.  Only  in  exceptional  cases  (e.g.  Poots*  (1965)  analysis  of 
the  compressible  stagnation  point  boundary  layer)  is  an  exact  compres¬ 
sible-incompressible  correlation  established,  while  compressibility  de¬ 
stroys  the  very  feature  (i.a. ,  the  independence  principle)  which  under¬ 
lies  many  of  the  available  incompressible  solutions. 
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XTI.  STABILITY  OF  THREE-DIMENSIONAL  LAMINAR  BOUNDARY  LAYERS 

The  theory  of  hydrodynamic  stability  of  three-dimensional  boundary 
layers  is  reviewed  by  Stuart  on  pp.  549-558  of  Chapter  IX  of  L. B.L.  An 
additional  important  survey  is  made  by  Brown  in  Vol.  II  of  B.L.F.C. 

The  starting  point  of  the  usual  theory  is  stated  by  Stuart  (in 
Gregory,  Stuart  and  Walker  (1955))  as  follows: 

"At  a  local  station  in  the  flow,  the  equations  of  sta¬ 
bility,  with  certain  approximations,  were  found  to  resemble 
formally  those  for  two-dimensional  flows;  the  relevant  mean- 
flow  velocity  is  the  component  in  the  direction  of  propaga¬ 
tion  of  the  disturbance  at  that  particular  station."* 

If  we  assume  that  the  spectrum  of  available  infinitesimal  distur¬ 
bances  is  isotropic  in  direction,  the  stability  theory  for  two-dimen¬ 
sional  flows  can  be  applied  to  the  projections  of  the  undisturbed  ve¬ 
locity  profile  on  various  planes  passing  through  the  local  wall  normal, 
to  find,  for  example,  the  mini. mum  critical  Reynolds  number  as  a  func¬ 
tion  of  the  angular  orientation  of  the  plane. 

The  most  prominent  result  of  such  a  study  (examples  are  given  by 
Brown  and  by  Gregory,  Stuart  and  Walker  (1955))  is  the  destabilizing 
influence  of  secondary  flew.  In  its  presence,  the  velocity  profile 
Inevitably  exhibits  an  inflexion  point  in  many  of  its  projections,  as 
shown  in  Fig.  7. 

Inflected  profiles  (see  for  example,  Schlichting  (1960)  p.  387) 
have  lower  minimum  critical  Reynolds  numbers  and  much  larger  areas  of 
instability  in  the  wave-number  versus  Reynolds -number  plane  than  do 
profiles  without  inflexion.  This  fact  partly  determines  the  appropri¬ 
ate  methods  of  analysis  of  the  stability  equation  (see  Brown  (1961)). 

An  infallible  rule  of  thumb  to  indicate  which  projection  of  the  ve¬ 
locity  profile  has  the  lowest  critical  Reynolds  number  does  not  seem 
to  exist,  but  from  examples  shown  by  Brown  it  seems  to  lie  close  to 
the  cross-flow  ($  =  90°)  profile. 


This  assumption  has  recently  been  criticized  by  Lilly  (1964) 
in  a  study  of  Eckraan  layer  stability  in  a  rotating  tank. 


/  » 


Fig.  7  —  Projections  of  the  velocity  profile 
on  planes  through  the  normal. 

A  projection  of  unique  importance  is  that  in  which  the  inflexion 
po„nt  coincides  with  the  crossover  point  (where  the  projected  u  is  ze¬ 
ro).  Such  a  profile  can  support  a  stationary  neutral  disturbance 
(phase  speed  =  0)  whereas  the  usual  neutral  disturbances  of  two- 
dimensional  boundary  layers  are  traveling  waves  with  C  >0.  (See, 
for  example.  Fig.  16.12,  p.  397  of  Schlichting  (I960)).  Such  sta¬ 
tionary  disturbances,  having  somehow  appeared  at  particular  locations 
on,  say,  a  swept  wing  or  a  rotating  disk  (Gregory,  Stuart  and  Walker 
(1955)),  manifest  their  continued  presence  by  a  street-like  pattern 
made  visible  by  china-clay  technique  even  substantially  upstream  of 
transition.  Stuart's  theory  successfully  predicts  the  "sweep  angle" 
of  these  streets  (which  lie  at  right  angles  to  the  projection  plane 
in  which  the  inflexion  and  crossover  points  coincide),  but  gives 
poor  estimates  of  the  observed  wavelengths.  Eichelbrenner  and  Michel 
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(1958)  have  made  observations  of  what  may  possibly  be  the  same  phenom¬ 
enon  occurring  on  a  flattened  ellipsoid.  Finally,  it  has  been  specu¬ 
lated  that  even  turbulent  layers  may  be  subject  to  this  sort  of  insta¬ 
bility,  evidence  for  which  is  the  regular  streets  of  Sargassum  or  foam 
occasionally  arrayed  over  large  areas  of  the  ocean  surface  (Faller, 
1964)  or  by  cloud  rows  in  the  atmosphere  (Faller,  1965). 
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I 

I  XIII.  TRANSITION  TO  TURBULENCE 


While  it  is  known  that  transition  to  turbulence  of  the  laminar 
boundary  layer  does  not  take  place  by  simple  amplification  of  two-di¬ 
mensional  disturbance  waves.  Smith  (1956)  was  able  to  find  a  useful 
description  in  terms  of  stability  theory  parameters  of  a  correlation 
of  smooth-wall,  low  free-stream  turbulence  transition  data  made  by 
Michel  (1952).  An  attempt  to  extend  the  Michel-Smith  transition  cri¬ 
terion  to  make  it  applicable  to  three-dimensional  boundary  layers  with 
weak  cross  flows  has  been  described  by  Eichelbrenner  and  Michel  (1958), 
and  tested  by  them  with  china-clay  observations  of  transition  on  a 
yawed  ellipsoid  of  revolution.  Because  the  comparison  seemed  lees  suc¬ 
cessful  than  one  might  have  expected,  particularly  along  the  flanks  of 
the  ellipsoid,  they  then  tried  a  second  transition  criterion,  that  of 
Owen  and  Randall  (1952  and  (1953),  which  appeared  to  be  related  to  the 
instability  of  the  secondary  flow. 

Owen  and  Randall's  criterion  is  derived  from  experiments  on  yawed 
wings,  and  does  not  pretend  to  any  great  generaltiy.  In  its  present 
state  of  development,  it  simply  associates  transition  with  a  single 
critical  value  of  a  Reynolds  number  based  on  boundary- layer  thickness 
and  the  maximum  value  of  the  cross-flow  velocity: 


Re 


tr 


max  v 


r  «  200  -  300 


where  the  range  of  values  presumably  comes  from  some  as  yet  unexplained 
dependence  on  profile  parameters. 

When  Eichelbrenner  and  Michel  plotted  points  of  v^^  6/u  «  300  on 
the  surface  of  their  ellipsoid  of  revolution,  they  found  much  better 
agreement  with  the  observed  transition  curve  at  overall  body  Reynolds 
numbers  U  L/v  ■  2  to  6  million  than  they  obtained  with  the  Michel- 
Smith  criterion.  Of  course,  only  the  latter  criterion  could  apply  on 
lines  of  symmetry  where  ■  0,  and  in  the  vicinity  of  such  lines  it 

might  be  presumed  to  be  the  preemptive  criterion. 

Finally,  Eichelbrenner  and  Michel  remark  that  their  observed  tran¬ 
sition  lines  frequently  exhibit  "tongues,"  the  number  of  which  appears 
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to  increase  with  increasing  Reynolds  number.  Similar  tongues  show  in 
china-clay  pictures  of  transition  on  swept  wings,  where  they  presumably 
develop  from  the  stationary,  swept  laminar  disturbances. 

While  the  "tongues”  of  Eichelbrenner  and  Michel  may  be  of  a  simi¬ 
lar  origin,  the  distances  between  them  are  not  negligible  compared  to 
radii  of  curvature  of  the  ellipsoid  itself,  so  that  one  of  the  basic 
approximations  of  the  standard  "plane-flow"  stability  theory  cannot  be 
satisfied. 

No  detailed  studies  of  the  mechanics  of  transition  of  boundary 
layers  with  cross  flow,  comparable  to  those  made  by  Schubauer,  Kleba- 
noff  and  others  for  the  case  of  two-dimensional  undisturbed  flows, 
seem  to  exist. 
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XIV.  THREE  -PI  MEKSI  ONA.L  TURBULENT  BOUNDARY  LAYERS 

On  exterior  aerodynamic  surfaces,  in  rotating  machinery,  and  in 
natural  flows  (rivers,  the  atmosphere,  etc.)  one  has  to  deal  with 
three-dimensional  boundary  layers  that  are  turbulent.  The  theory  of 
such  boundary  layers  is  only  very  slightly  developed,  and  perhaps  be¬ 
cause  of  the  formidable  aspect  of  even  the  laminar-flow  differential 
equations  there  has  been  almost  no  work  along  deductive  lines,  utiliz¬ 
ing  ad  hoc  models  of  eddy  diffusivity  or  mixing  length  to  provide  the 
missing  link  between  Reynolds  stresses  and  mean-velocity  derivatives. 
Exceptions  can  be  found  in  the  meteorological  literature,  but  they 
typically  involve  special  assumptions  (negligibility  of  boundary- layer 
growth  rate  or  of  horizontal  mean-flow  convection  which  set  them  fairly 
far  apart  from  the  main-stream  of  boundary- layer  theory. 

Most  studies  of  turbulent  three-dimensional  boundary  layers  have 
thus  far  taken  an  inductive  and  empirical  approach,  seeking  out  ways 
by  which  measured  mean- velocity  profiles  can  be  represented  parametri¬ 
cally,  in  hopes  that  momentum-integral  methods  can  be  devised  to  pre¬ 
dict  the  evaluation  of  profile  parameters  for  given  bodies,  surface 
conditions  and  external-flow  velocities. 


TIME -AVERAGED  BOUNDARY -LAYER  EQUATIONS 


The  boundary- layer  approximations  to  the  time-averaged  Navier- 
Stokes  equations,  for  the  case  of  constant  fluid  density  and  viscos¬ 
ity,  and  for  a  steady  but  three-dimensional  tirae-averaged  flow,  may 
be  derived  by  the  following  sequence  of  assumptions. 

The  boundary  layer  is  thin  by  comparison  with  any  radius  of  cur¬ 
vature  of  the  body  surface,  so  that  we  may  still  assume  h^  **  1, 
dh^/d£  >=>  dh^/SC  =  0  in  the  region  of  interest.  With  these  simplifi¬ 
cations  the  instantaneous  momentum  and  continuity  equations  can  be 
combined  to  give 
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3u  ,  1  S  ,2.  ,19  .  x  ,  9  ,  v./2  2.  ,  _ 

a?  +  h^H^u)  +  h;Si  (uv)  +  5j  <“»)  +  <«  -  »  )»2  +  2“v>li 


.  ,  (  9  ?\  /'•-  other  viscous  terms  which  \ 

—  Te(A  +  ^-  —  tuR  ]  -  2u)  v  =  v  i  — ^  +  are  linear  in  the  velocity/ 
1  o C  components 


fi  +  ^  li  (vu>  +  ^  In  <v'2)  +  §c <TO)  +  (v2 '  “V  +  2vuk2 


^  -  ~  u>2R2)  +  2in.jU  =  u  0-^  +  other  linear  vi 
2  oC 


viscous  terms 


+  i_  (wu)  +  JL  (w)  +  —  (w2) 

9t  h1  9§  '■wu;  +  h2  9T,  KVV}  9C  ^  ; 


.  0  ( .  ,  p  1  2  2\  ,  0/  .  f3w  ,  other  linear  ^ 

+  dcVA+p2U,Ri  +  2(t0lV  "  ®2U>  SV  \^~2+  Viscous  terms/ 


In  these  equations 


1  ahi 

‘l  '  hjh2  oil  • 


and  h 


t  8h2 
2  "  hxh2  9| 


are  the  geodesic  curvatures  of  the  surface  coordinate  curves. 

Into  these  equations  we  now  introduce  the  decomposition  of  tur¬ 
bulent  variables  into  mean  and  fluctuating  parts,  and  then  average 
over  time.  For  example,  we  write  u  =  u  +  u7,  with  the  understanding 
that  u  =  u,  u7  =  0.  Then 


-2  -2  0-  /  ,2  -2  ,  ,2 

U  »  U  +  2UU  +  U  =  U  +  U 
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We  assume  that  du/dt  =  0.  To  save  space,  ve  simply  Imagine  this  oper¬ 
ation  to  be  carried  out. 

We  now  assume  that  the  thinness  of  the  layer  implies  that  if  1/h^ 

3/5§(  )  and  l/l^  d/dT|(  )  are  taken  to  be  of  order  unity,  then  d/d£(  ) 

i3  an  order  of  magnitude  larger  (say  0(L/6)  ,  where  6/L  «  1).  The 

q.'js&tixy  in  (~)  may  be  any  averaged  property  of  the  velocity  field, 

like  u  or  (u^).  For  the  moment  we  use  this  assumption  only  to  estab- 

2-  2 

lish  the  fact  that  the  mean-flow  viscous  terms  (e.g.,  u3  u/d£  ),  re¬ 
tain  the  same  form  as  in  a  laminar  boundary  layer,  and  to  establish, 
from  the  time -averaged  continuity  equation,  that  if  u  and  v  are  taken 
to  be  of  order  unity,  then  w  is  of  order  6/L  «  1.  With  these  approx¬ 
imations  the  time-averaged  boundary-layer  equations  become 


^  Is  ("2>  + 12 1 if  ("*>  +  k  <“”>  +  (a2  -  ;2>k2 

+  2“S  +  (A  +  p  -  \  “2r2)  -  V 


^  I?  +  If  (’)Z  +  fc (iS)  +  (’2  ‘  “2)si 


+  2vuh.  +  — 
2  h„ 


(A  +  i .  I  A2)  + 
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and 


Sc  (A  +  P  ■  1  ,1,2r2)  +  2(V  '  U,2G)  = 


a_ 


(v'u^  +  JL  iL  (w  »v'\  +^- 
^  u  '  +  h2  ST]  V  }  dC 


To  make  further  simplifications  we  must  have  additional  experi¬ 
mental  knowledge,  regarding  the  typical  orders  of  magnitude  of  the 
kinematic  Reynolds  stresses  (e.g. ,  u7  ,  uv7)  and  their  derivatives. 
This  information  is  not  directly  available  for  three-dimensional  tur¬ 
bulent  boundary  layers,  so  we  shall  have  to  employ  a  liberal  amount 
of  guesswork.  We  should  note  with  particular  care  that  while  v7  and 

O 

u7  are  probably  roughly  independent  of  our  choice  of  surface  coordi¬ 
nate  system  (in  a  two-dimensional  boundary  layer,  streamwise  fluctua- 

-2 

tion  and  cross-stream  fluctuation  are  roughly  equally  energetic),  v 
-2 

and  u  do  depend  on  this  choice.  We  shall  proceed  in  intrinsic  coor¬ 
dinates,  and  shall  assume  that  all  kinematic  Reynolds  stress  compo¬ 
nents  are  of  the  same  order  of  magnitude,  which  vc  estimate  (from  Kle- 
banoff’s  data  for  a  zero  pressure  gradient  two-dimensional  boundary 
layer)  to  be  no  more  than  1  percent  of  the  local  U  .  Thus  we  shall 
drop  1/h^  d/B§  (tPv7)  and  l/h^  9/ST)  (v7)^  compared  to  d/d£  (v7 w7)  in 
the  T|-momentum  equation,  and  make  corresponding  simplifications  in  the 
other  equations.  The  ^-momentum  equation  can  be  integrated  over  Q  to 
give 


PCSill.C)  =  A  +  ^  ~  2  w2r2  “  P0<M,«>)  -  w'2  +  2  J  (uyi  -  (OjV)  dQ 

o 

With  the  assumption  that  any  coordinate  spin  is  moderate,  so  that 
ojU6  is  at  most  of  order  6/L  «  1,  we  can  drop  the  integral.  Then,  when 
»/a§  and  BP  /BT)  are  evaluated  from  the  inviscid  momentum  equation,  we 

°  °  2  _2  2 

shall  be  comparing  3/35  (U  )  or  u  with  3/35  w7  ,  and  we  assume  that 
the  fluctuation  term  can  be  dropped  in  that  comparison.  Substituting 


for  dPQ/d§  and  dPo/dT}  in  the  §-  and  ^-momentum  equations,  and  using 
the  continuity  equation  again,  we  find  the  forms 


u  du  U  3U  ,  v  tlu  •  5u  -2  .  —  . 

ri7TT?  +  rTr'i'w  V  *  V  h  +  UVH  -  2uj.v 
hL  d|  hx  de,  h2  dC  5C  2  13 


U  “  I?  "  (u#2  -  v/2)%  -  2u  #v  V 

dC2  dC  2  1 


u  dv  ,  v  dv  ,  -  du  .  .,,2  -2n  .  —  .  „  - 

a?  +  Si +  “  ac  (  +  2l”3u 


u  (v^S7)  -  (v'2  -  u'2)Kj  -  2u'v  'k2 


We  reserve  judgment  on  the  importance  of  the  last  two  terms  in 
each  equation  until  we  have  integrated  these  equations  over  Q< 

MOMENTUM-INTEGRAL  EQUATIONS 


We  note  the  boundary  conditions  at  £  **  0, 


u  =  v  =  0 


(no  slip) 


w  =  w 


(given  blowing  velocity) 


U  W  “  V  w 


(for  a  smooth  wall) 


du  — t — ; 
u  ac  -  u  w 


2  Cf, 


and 


dv  ~r~. 

v  K  "  v  w 


£c. 


(at  a  rough  wall) 


AC  C  = 


u  =  U,  v  »  0 


ou  5v 

SC  =  sc 


v  'w !  =>  0 


We  define  the  integral  thickness  as 


s  f(l -{})«• 


CO  - 

5  -J 


l0n  B  J  5  (x  -  5)  «> 


Lei2 =  r  §  g  - 1) 


-®  -l 


=  -f§«  - 


Next  we  add  (u  -  U)  times  the  time-averaged  continuity  equation  to  (10) 
and  v  times  the  continuity  equation  to  Eq.  (11) .  This  gives 

/u  -  iA  au  i  a  i  ^  v  au  ,  a  r-,- 

5?  +  iT[  5?  [u(u  •  U)1  +  iq  si)  [v(“  ‘  “>  1  +  ^  sii +  ac  [u(u  - 


-  ....  .  >2- 
+  [(u  -  U)u  -  v  ]h  +  [(u  -  U)v  +  uv]k  -  2u)_v  =  v  — ~ 

J  ac 


|p  (u^7)  -  (u/2  -  V,2)K2  -  ?,U  /V,K1 


•■-is 


I xl  tTi> ?  UfJ»A>* i-'AV<  4  VWVj •NtU'tVLfc . .> -'JA'-'i  *\ '.,WD 1A«* U VAS%j4i  (f*/> *  »li  ?UZ  <* .  kIU^4UHiUi)'Uv(U«^>v..  ,,  t iw’.'W  AUrl  J  V,  1  fl » WiUtfi S. «Tlii  L  ilMllt 
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and 


^ tc (55) +'S2' li  (’>2  +  tc (;s> 


(u2  -  -2  •  -2 


u  4-  v  )x^  +  2uvh2 


I  ' 


+  2w  (u  -  U)  *  y  (v  lv ')  -  (v'2  -  u/2>*  -  2u  Vv  7h 

3  £^2  o£  12 


All  terms  will  now  give  finite  values  when  integrated  from  £  «  0  to 

C  **  o.  (In  these  integrations  U  and  h.. ,  h  are  taken  independent  of 

2  12 

£. )  After  division  through  by  U  and  introduction  of  the  symbols 


x  .  I  v  =  a 

X  L’  y  L 


(and  other  symbols  such  as  m,  n,  X,  h,  as  used  above),  ve  obtain 


1  *911  1  3012  fl  1 

TT  * +  r2  IT  =  -T  -  h^  W2fln  +  V  +  K(en  -  e22  -  Fn):> 


2X 

h2y 


20  L  /2u)»  4-  ILv  w 

F12  +  U  012  +  V  U  )  L62  +  U 


2u),  +  0- 


and 


1  a021  1  "22  f2  2 

+  22  „  — i  .  -2-  [(m  +  H)e -  xF  10] 

hj  dx  h2  dy  2  h^x  LV  21  12J 


tuy  ^X(022  *  6li  *  Fll)  +  n6l^ 


20- L 


+  0„ 


+  _1I  e  f  3  " 

+  U  022  "  V  U  )  L61 


-75- 


Note:  Except  for  the  presence  of  the  blowing  velocity  and  the  "fluc¬ 
tuation  integrals"  and  these  equations  agree  exactly  with  the 

laminar  flow  Eqs.  (5)  and  (6).  The  apparent  differences  in  their  re¬ 
spective  right-hand  numbers  came  from  the  fact  that  in  our  laminar-flow 
treatment  we  have  measured  6's  and  9's  in  units  of  z,  where 
a  -^UL/ux(£/L),  while  in  the  present  treatment  they  are  measured  sim¬ 
ply  in  units  of  Qfh.  This  means  that 

lx  911,21(z)  =  3x  911,21  (l)  *  2x~  011,21  (l) 

3y  812,22(z)  =  by  912,2 2.  (l)  +  2y  9 12,22  (l) 

Similarly  we  can  show  that 


f"(z  -  0) 
eu<*> 


In  their  present  fora  (with  6's  and  0‘s  measured  in  units  of  £/L)  and 

with  the  suction- velocity  term  w  /U,  the  equations  arc  easily  special- 

s 

ized  back  to  their  familiar  forms  for  two-dimensional  boundary  layers 
(h^  “  1,  h  “  X  3  fij  m  u?2  w  0).  Then  the  streamwisa  momentum- integral 
equation  becomes 


+  v!  + 


w 

8 


u 


ELI  HI  NATION  OF  THE  FLUCTUATION  INTEGRALS 

If  we  can  make  order-of-magnitude  estimates  from  data  on  two- 
dimensional  boundary  layers,  it  appears  as  though  F^  may  not  exceed 
2  percent  of  9^,  and  hence  it  may  be  dropped  wherever  the  two  quanti¬ 
ties  appear  as  a  sura  or  difference.  (This  relies  on  8^  +  8  ^  being 
not  much  smaller  than  9^,  an  assumption  which  could  only  be  violat  ed 
in  case  of  extraordinarily  strong  cross  flows.) 
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We  cannot  so  decisively  dispose  of  by  comparison  with  9^  or 

s^nce  it  is  hard  to  decide  on  typical  orders  of  magnitude  for 
either  quantity.  Since  further  development  of  the  theory  would  be 
nearly  impossible  in  a  case  in  which  F^  was  important,  we  shall  ar¬ 
bitrarily  ignore  it.  The  equations  which  we  shall  discuss  further  are 
then  cast  entirely  in  terms  of  integral  quantities  which  can  be  de¬ 
rived  from  assumed  profiles  of  u  and  v  and  the  friction  coefficients 

C,  and  Cf  . 

1  £2 

COUNT  OF  UNKNOWNS  AND  EQUATIONS 

If  we  accept  the  approximations  whereby  "fluctuation  integrals" 
such  as  F^  and  F^  are  neglected,  we  are  left  with  two  differential 
equations  for  the  determination  of  seven  unknowns  (five  independent 
integral  thicknesses  and  two  skin  friction  coefficients). 

For  subtle  reasons  which  could  hardly  be  anticipated  before  an 
examination  of  experimental  u  and  v  profiles  (see  Rotta  (1962)  PP- 
70-73  and  pp-  172-181),  the  compatibility  conditions  at  the  wall  do 
not  provide  useful  additional  relations  between  the  6's,  8's,  and  say, 
the  given  pressure  gradient. 

Additional  weighted  momentum-integral  equations  can  of  course  be 
formed,  but  we  can  see  at  a  glance  that  these  will  involve  important 
integral  thicknesses  in  which  the  Reynolds  stresses  appear  in  the  in¬ 
tegrands.  For  example,  the  u-weighted  ^-momentum  equation  will  involve 
the  thickness 


Jo  ft  ^ dc 

which  describes  part  of  the  conversion  of  kinetic  energy  from  the 
"mean-flow  budget"  to  the  "fluctuation  budget."  Since  the  Reynolds 
stresses  are  harder  to  measure,  with  acceptable  accuracy,  than  the  u 
and  v  profiles  from  which  the  6's  and  0's  are  computed,  there  does  not 
yet  exist  an  extensive  empirical  base  for  work  with  the  weighted  inte¬ 
gral  equations.  If  we  abandon  them,  we  are  left  to  devise  five  addi¬ 
tional  relations  between  the  6's,  8's,  Cf's  and  our  given  parameters, 
by  direct  examination  of  experimental  data. 
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EMPTRICAL  GENERALIZATIONS  ABOUT  TWO-DIHENSI  ONA.L  TURBULENT  BOUND  ARY  - 
LAYER  PROFILES 


In  warty  cases,  especially  those  with  weak  cross  flows,  we  expect 
the  streamwise  velocity  profile  to  be  quite  like  those  which  have  been 
so  extensively  studied  in  flows  possessing  a  carefully  cultured  two- 
dimensionality  of  the  time-averaged  state.  The  properties  of  such 
flows  have  been  reviewed  in  detail  recently  by  Rotta  (1962)  and  Schu- 
bauer  and  Tchen  (1959).  We  are  mostly  interested  in  the  empirical 
generalizations  for  u  which  are  reviewed  by  Rotta  on  pp.  156-166.  In 
particular,  we  select  for  further  discussion  the  doubly-infinite  family 
of  u  profiles  given  by  the  wall-wake  model  of  Coles  (1956). 


The  Wall-Wake  Model  of  u  Profiles 

According  to  this  model,  which  though  frankly  empirical  and  ap¬ 
proximate,  is  capable  of  fitting  experimental  u  profiles  with  remark¬ 
able  accuracy  throughout  a  wide  range  of  both  favorable  and  adverse 
pressure  gradients,  the  profile  is  constructed  by  a  linear  combination 
of  two  universal  functions.  One  of  these,  called  the  wall  function  or 
the  law  of  the  wall,  relates  u  only  to  the  local  friction  velocity. 

To,  the  kinematic  viscosity  v,  and  the  normal  distance  C>  pro- 
T  W 

vided  that  the  wall  is  smooth.  For  nonzero  skin  friction,  this  law 
dominates  the  u  profile  over  some  finite  range  of  £,  and  within  this 
region  the  local  pressure  gradient  has  no  direct  effect  (within  the 
broad  range  of  conditions  on  which  the  law  is  based).  By  dimensional 
analysis,  this  law  can  be  written 


u 

u 

T 


When  the  wall  is  rough,  and  if  the  roughness  elements  possess  geometri¬ 
cal  similarity  in  their  shape  and  spacing  so  that  they  can  be  charac¬ 
terized  by  a  roughness  height  k^,  the  law  may  be  generalized  to 
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u 

u 

T 


In  some  cases,  particularly  those  involving  favorable  pressure 
gradients  or  large  Reynolds  numbers,  the  wall  function  may  suffice  to 
describe  u  over  a  sufficient  range  of  ut  Q/v  so  that,  for  u^  £/u  >  50, 
it  assumes  the  famous  logarithmic  form 


with  K  “  0.41  and,  for  smooth  walls,  C  a  5.0,  as  found  from  experiment. 
When  thi3  is  attained  for  a  rough  wall,  the  effect  of  roughness  is  en¬ 
tirely  in  the  "constant1*  C: 


£  in  (~u“)  +  0  (“V1,  roughness  shape)  (ut  £  >  50,  Q>  kr) 


The  dependence  of  C  on  u^k^/v  is  independent  of  shape,  but  for 
Nikuradse's  sand  grains,  u  k  /v  £5  gave  C  =  5.0,  or  smooth  flow, 
while  for  u^/v  £  70,  fully  rough  flow  with  C  ■  (shape  only) 

(1/K)  Xn  (u^k^/v),  so  that  u/u^  =  (1/K)  Xn  (C/k^)  +  C^  (fully  rough). 

The  direct  influence  of  local  pressure  gradient  and  boundary- 
layer  thickness  enter  u  through  the  second,  or  wake,  component  of  the 
profile.  The  form  of  this  was  discovered  by  Coles  (1956)  by  a  proce¬ 
dure  indicated  in  Fig.  8.  Only  profiles  for  which  the  logarithmic 
portion  of  the  wall  law  was  empirically  well  established  could  be 
used  initially.  The  wake  component  was  determined  by  subtracting  the 
extrapolated  (logarithmic)  wall  law  from  the  measured  u,  and  a  defi¬ 
nite  procedure  for  defining  the  boundary  layer  thickness  6  was  estab¬ 
lished  (as  seen  in  Figc  8).  Ideally  u  as  well  as  u  should  be  direct- 

T 

ly  measured  in  the  experiment. 

The  experimental  residual,  u  *  u^f,  was  discovered  by  Coles  to 
have  a  universal  shape  in  boundary  layers  under  a  free  stream  of  suf¬ 
ficiently  low  (<  0.1  percent)  turbulence  level.  With  6  defined  by 


Midpoint- 


Measured  points- 


Tangent  to  data, 
parallel  to  I“1 

II 


Line  jj  to  I  &  II 
half-way  between 


jrithmic  portion  o 2  the 
!  law  of  the  wall 


j  /—This  value  of  £  =  —  by 
\S  assumption 


£  (actual  measured  distance  from  wall, 
plotted  on  a  logarithmic  scale) 

Fig*  8  --  Law  of  the  wake. 

Fig.  8,  and  with  a  universal  wake  function  W(£/5)  normalized  to  run 
from  0  to  2  as  C  goes  from  0  to  6,  Coles'  complete  wall-wake  modal  for 
u  is  given  by 


,u  Q  uk  \ 

■  £i  (Hr-  -VO  +  £  “  ©•  0  £  y  s  6 


In  this  the  dimensionless  parameter  n(|)  assumes  values  determined  by 

the  local  pressure  gradient  and  prior  boundary- layer  history,  and  K 

is  again  the  von  ICarman  constant.  We  can  re-express  it  in  terms  of 

the  given  U  (instead  of  the  unknown  u  )  as 

T 

l  - V5 (t) 


We  see  that  this  forms  a  doubly-infinite  family  of  profile  shapes, 
governed  by  the  three  dependent-variable  parameters  6/L,  C^,  and  tt  and 
the  two  given  parameters  UL/u  and  k^/L.  By  setting  £  =  6  we  get  one 
relation  between  these  parameters,  namely 


This  difficult  implicit  relationship  is  untangled  graphically  by  Rotta 
(1962)  on  pp.  168  and  169  (for  C  =  5*2). 

The  integrals  for  61  and  0^  can  now  be  carried  out  by  noting 
that  K(0)  d0  =  1,  and  ignoring  deviations  of  f^  from  its  asymptotic 
logarithmic  form,  so  that  £  Sf^/3£  d(£/6)  «  1.  We  get 

6i  -  6V5(4ji) 

and 

9u  -  8  {^f  (Hr1)  -  +  o'"  +  *>"2> 

where,  according  to  Coles,  ot  =  1.600  and  &  =  0.761. 

The  "Missing  Relation” 

When  6^  and  9^  are  eliminated  in  favor  of  6,  C^,  and  rr,  the  mo¬ 
mentum-integral  equation  plus  the  free-stream  relation  (boxed  above) 
give  us  two  equations  for  three  unknowns. 

Neither  of  these  equations  postulates  any  direct  local  connection 
between  pressure  gradient  and  profile  shape,  and  we  might  hope  to  find 
something  like  this  empirically,  perhaps  in  the  form  of  a  correlation 

between  rr  and  6.(dp/d§)/T  •  Each  of  these  quantities  is  constant 
i  o 

(Coles  1956,  Clauser  1956)  in  so-called  equilibrium  or  self-preserving 
boundary  layers,  so  that  a  few  points  on  a  correlation  curve  are  readily 
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available.  Of  course,  it  is  naive  to  hope  that  a  local  correlation 
of  this  type  might  hold  very  generally  in  boundary  layers  in  which  tt 
changes  radically  with  §,  but  even  a  rough  correlation  would  probably 
lead  to  results  as  good  as  those  obtainable  with  other  popular  meth¬ 
ods — e.g.,  that  of  von  Doanhoff  and  Tetervin  (see  Schubauer  and  Tchen, 
1959)*  Good  accuracy  would  be  expected  in  favorable  pressure  gradients 
(where  tt  is  small,  rang: ng  between  zero  in  converging  channel  flows  to 
0.55  with  zero  pressure  gradient),  and  poorer  accuracy  would  be  ex¬ 
pected  as  dp/dx  becomes  increasingly  positive  and  separation  is  ap¬ 
proached.  In  the  latter  case  the  detailed  history  of  boundary- layer 
growth  is  bound  to  play  an  increasingly  important  role,  even  more  so 
than  in  laminar  layers.  In  the  extreme,  Clauser  (1956)  discovered 
that  the  profile  in  strong  adverse  pressure  gradients  was  so  sensitive 
to  upstream  history  that  it  could  not  be  effectively  manipulated  by 
changing  the  local  dp/d§.  This  phenomenon  is  occasionally  referred 
to  as  instability,  in  the  sense  that  a  small  local  change  in  dp/d§  (or 
in  wall  roughness)  may  produce  a  downstream  effect  which  grows,  rather 
than  diminishes,  with  increasing  downstream  distance. 

Finally,  we  note  with  Rotta  (1962),  p.  199,  that  the  "missing 
relation,"  as  we  have  called  it,  must  be  an  independent  empirical  dis¬ 
covery,  and  that  it  cannot  be  obtained  by  any  manipulat?„on  of  the 
equations  (momentum  and  continuity)  and  profile  information  (d(C/6, 

6/L,  C^,  tt))  we  already  possess.  For  example,  even  though  Coles  has 
shown  that  the  wall-wake  model  can  be  applied  to  the  momentum  equa¬ 
tion,  and  accurate  values  of  Reynolds  stresses  computed,  the  result¬ 
ing  stress  profile  cannot  be  inserted  into  (say)  the  energy  integral 
equation  to  obtain  our  missing  relation.  What  results  from  such  a 
procedure  is  an  identity  ( 1  =*  1)  and  not  an  independent  equation. 

EMPIRICAL  MODELS  FOR  THREE-DIMENSIONAL  TURBULENT  BOUNDARY-LAYER 
PROFILES 

We  found  in  our  review  of  the  two-dimensional  turbulent  boundary 
layer  that  there  exists  a  profile  model,  the  wall-wake  model,  which 
accurately  fits  a  very  wide  variety  of  profiles  by  use  of  two  univer¬ 
sal  functions  and  three  parameters,  even  though  we  do  not  yet  possess 
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an  effective  way  of  predicting  these  parameters.  We  now  wish  to  see 
whether  in  the  three-dimensional  case  we  can  go  equally  far,  discover¬ 
ing  widely  applicable  models  for  both  the  streamwise  and  cross-flow 
velocity  profiles.  As  might  be  expected  from  the  comparatively  small 
amount  of  study  yet  given  to  three-dimensional  boundary  layers,  the 
answer  seems  to  be  "not  yet."  Two  major  proposals,  both  of  which  seem 
to  have  some  experimental  support,  have  been  made;  they  agree  in  some 
cases  and  disagree  in  others. 

Generalization  of  the  Wall-Wake  Model 

In  his  1956  article  Coles  suggested  a  straightforward  vector  gen¬ 
eralization  of  the  wall-wake  model  into 


-  -  f 
u  =  u  f,  l - , 

T  1  \  V 


x  k  v  u 

+  n  -£ 
V  /  -  K 


where  f^  and  W  are  the  same  universal  functions  as  in  the  two-dimen¬ 
sional  case,  and  n  (C,T|)  has  become  a  tensor  parameter,  which  operates 
on  the  friction  velocity  to  beget  another  vector  tangential  to  the 
wall.  The  friction  velocity  vector  has  the  direction  of  the  skin  fric- 


Experimental  testing  of  this  model  requires  accurate  measurements 
of  both  magnitude  and  direction  of  velocity  and  skin  friction,  and  it 
appears  that  instrumentation  employed  so  far  is  somewhat  cruder  than 
that  used  in  two-dimensional  flow  studies.  Be  that  as  it  nay,  the 
test  would  in  principle  consist  of  two  steps:  (a)  a  check  of  the  par¬ 
allelism  of  the  wake-component  vectors,  as  visualized  in  the  polar 
plot  of  Fig.  9;  (b)  a  plot  of  the  wake- component  magnitude  versus  £. 

If  the  data  pass  the  first  test  (parallelism  of  wake-component  vec¬ 
tors)  we  then  see  whether  the  amplitude  of  the  wake  component  can  be 
described  by  a  universal  profile  function.  Finally,  if  this  wake  pro¬ 
file  shows  the  symmetry  of  W(£/6),  6  is  easily  determined,  and  so  is 

•■4 

the  amplitude  of  the  vector  tt  u  . 

Since  1956,  various  experimenters  have  compared  data  with  the  gen¬ 
eralized  wall-wake  model;  while  none  of  their  comparisons  quite  live  up 
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Many  polar  plots  were  found  to  be  nearly  triangles,  with  base  given 
by  the  inviscid  velocity  vector  U,  one  side  aligned  along  the  skin 
friction  direction,  and  the  third  making  an  angle  with  U  which  appeared 
in  many  cases  to  be  simply  related  to  the  total  angle  through  which  the 
inviscid  streamlines  had  turned  after  a  point  in  which  the  flow  had 
originally  been  two-dimensional.  The  triangle  law  is  thus  a  simple 
statement  of  the  dependence  of  v  on  u,  without  direct  reference  to  Q. 

It  is  a  two-part  statement.  In  Part  I,  which  is  clearly  related  to 
Coles'  wall -dominated  region. 


v  =  eu 


In  Part  II,  the  outer  region, 

v  =»  A(U  -  u) 

The  corner  at  which  I  and  II  join,  was  estimated  originally  by 
Johnston  to  fall  at  u^/u  ~  16,  but  much  higher  estimates  (~  150)  are 
given  by  Hornung  and  Joubert  (1963). 

Conflict  Between  the  Two  Models 

While  certain  sets  of  data  (for  example,  that  of  Kuethe  et  al. 
(1949),  which  Coles  first  examined)  seem  to  fit  both  models  fairly 
well,  it  can  readily  be  seen  that  the  two  models  are  in  fairly  di¬ 
rect  conflict.  For  example,  we  can  show  that  if  Coles'  wake-compo¬ 
nent  vectors  are  indeed  parallel,  and  if  Johnston's  triangle  rule  is 
rigorously  applied,  the  wake-component  magnitudes  cannot  be  described 
by  an  S-shaped  function  like  w(|/6).  Under  the  assumed  conditions, 
the  angles  a  and  P  in  the  polar  plot  (Fig.  11)  are  independent  of 
C«  Thus  we  get 


wake 
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sin  or 
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Fig.  11  —  Hypothetical  projected  velocity  profile. 


If  we  assume  that  the  wall  law  is  logarithmic  for  £  >  and  define 

W  -  2u  .  (0/u  ,  (A),  we  find  that 

wake  wake 

w  =  2  in  (£-)/in  ft-) 

Curiously,  this  gives  a  fair  fit  to  Hornung  and  Joubert's  wake  func¬ 
tion  data,  as  shown  in  their  Fig-  7,  with  the  reasonable  choices  of 
6/£^  =  20  or  30 l  Coles  has  suggested  that  the  generalized  wall-wake 
model  is  most  likely  to  succeed  in  cases  in  which  the  divergence  of 

Reynolds  stress  is  reasonably  large  compared  with  the  lateral  compo- 

2 

nent  of  inertial  force,  i.e.,  when  d-r/dC  C  pu  k. ,  and  he  estimates 

2  1 
pu  ^/(dr/dQ  to  be  about  250  h^6.  If  the  radius  of  curvature  of 

surface  streamlines  is  much  shorter  than  a  few  hundred  boundary- 

layer  thicknesses,  the  outer  part  of  the  boundary  layer  accelerates 

as  a  more  or  less  inviscid,  although  rotational,  flow.  The  latter 

situation  has  been  assumed  by  various  authors  to  be  the  one  for  which 

the  triangle  model  will  be  most  successful. 

SUMMARY 

One  could  go  on  to  point  out  sample  flows  for  which  Johnston's 
triangle  rule  cannot  conceivably  apply,  but  perhaps  it  is  clear  from 
what  has  been  already  said  that  the  prognosis  for  successful  computa- 
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tions  of  three-dimensional  turbulent  boundary  layers  is  quite  poor  at 
present.  There  are  undoubtedly  problems  involving  favorable  pressure 
gradients  and  gentle  streamline  curvatures  in  which  success  nearly 
comparable  to  that  available  in  two-dimensional  problems  is  possible. 
For  example,  Eichelbrenner 's  (1963)  extensive  computations  for  ring 
wings,  and  those  of  Vaglio-Laurin  (1959)  for  reentry  bodies,  use  very 
simple  profiles  and  skin-friction  laws.  Turbulence  complicates  the 
problem  in  a  fundamental  way,  but  it  also  probably  suppresses  anoma¬ 
lous  cross  flows  by  its  vigorous  momentum  exchanges,  and  qualitatively 
am  be  expected  to  have  a  beneficial  effect  upon  the  prevention  or  de¬ 
lay  of  separation. 

It  seems  that  detailed  experimental  work  on  three-dimensional 
turbulent  boundary  layers  is  not  only  needed  in  its  own  right,  but 
may  also  somenow  help  us  to  sift  hypotheses  or  generate  fresh  ideas 
about  the  behavijr  of  two-dimensional  flows.  The  problems  of  instru¬ 
mentation  for  such  studies  are  quite  acute,  and  it  is  even  more  diffi¬ 
cult  to  attain  high  Reynolds  numbers  and  desirably  thick  boundary  lay¬ 
ers  than  in  studies  of  two-dimensional  flows. 
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XV.  SUGGESTIONS  FOR  FUTURE  WORK 

These  suggestions  do  not  reflect  a  review  in  depth  of  all  important 
and  outstanding  problems  of  the  field,  but  represent  a  collection  of 
items  which  appeared  to  be  most  intererting  and  perhaps  significant 
within  the  scope  of  this  study. 

THEORETICAL  TOPICS 

1.  Davey's  "strong  saddle  points  of  attachment,"  and  their  place 
in  a  complete  flow  field.  These  solutions  exhibit  a  kind  of  "reverse 
flow  without  separation,"  or  "harmless  separation."  When  situated  on 
a  leading  edge  between  two  adjacent  nodes  of  attachment,  the  solution 
for  the  singular  region  implies  a  fascinating  flow  structure  in  which 
the  saddle  point  of  the  inviscid  flow  field  is  transformed  into  some¬ 
thing  like  an  unusually  simply-described  reattachment  node  in  the  sur¬ 
face  flow.  Whether  this  really  happens  in  nature  is  apparently  an  open 
question,  and  a  convincing  theoretical  clarification  would  probably  be 
heralded  as  a  tour  de  force  in  boundary- layer  theory. 

2.  Further  development  of  the  three-diraensionalization  of  Smith 
and  Clutter's  numerical  procedure.  What  lias  been  given  above  is  only 
a  very  preliminary  sketch.  Seme  careful  analysis  might  be  worthwhile 
of  finite-difference  procedures  in  the  surfaces  of  constant  z,  of  nu¬ 
merical  stability  and  error  control,  and  some  computations  for  compari¬ 
son  with  the  method  of  Raetz. 

3.  The  method  of  characteristics  for  solving  momentum-integral 
equations.  Mager  has  brushed  off  as  impractical  Timman's  original  sug¬ 
gestion  that  the  momentum-integral  equations  be  integrated  along  their 
real  characteristics  because  the  two  characteristic  directions  include 
too  small  an  angle  in  the  case  of  weak  cross  flows.  However  true  this 
may  be,  it  does  not  seem  proper  to  close  the  subject  unless  one  wishes 
from  the  start  to  abandon  momentum- integral  methods  in  problems  with 
large  and  interesting  cross  flows. 

Even  for  the  weak  cross-flow  case  it  seems  to  me  that  the  method 
of  characteristics  might  be  employed  in  the  formulation  of  an  approxi¬ 
mate  integration  method  which  might  be  more  correct  and  efficient  than 
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that  currently  employed.  Extending  the  analysis  to  cover  the  method 
employing  both  unweighted  and  weighted  integral  equations  would  be  a 
challenging  investigation. 

EXPERIMENTAL  TOPICS 

1.  Mapping  of  skin  friction  field  in  laminar  flow.  Serious  ef¬ 
forts  have  been  made,  particularly  in  France  and  England,  to  observe 
skin  friction  trajectories  by  oil  and  lampblack,  or  by  dye  emission 
techniques.  The  results  have  been  very  helpful  and,  in  some  cases, 
almost  definitive.  It  would  be  a  major  achievement  if  these  or  other 
techniques  could  be  further  perfected  and  employed  to  generate  a  de¬ 
finitive  portfolio  of  experimental  skin- friction  fields,  including 

in  particular  some  which  could  provide  quantitative  checks  on  theory, 
and  some  for  pedagogical  illustration  of  the  richness  of  qualitative 
possibilities. 

2.  Yawing  of  turbulent  wakes.  A  further  understanding  of  the 
response  of  the  wake  component  of  turbulent  shear  flows  to  gradual  and 
to  sudden  yawing  might  contribute  significantly  to  the  construction  of 
better  models  of  the  three-dimensional  turbulent  boundary  layer.  A 
flat-plate  wake  might  be  studied  as  the  wind  tunnel  downstream  of  the 
plate  executes  a  bend  around  an  axis  normal  to  the  "plane"  of;  ti:e 
plate.  This  would  provide  an  opportunity  for  the  refinement  of  direc¬ 
tion-sensing  instrumentation  in  an  environment  that  is  less  cramped 
than  the  boundary  layer. 

3.  Synoptic  exploration  of  three-dimensional  turbulent  boundary 
layers.  Careful  synoptic  measurements  of  mean-velocity  profiles,  pres¬ 
sure  distributions  and  surface  stresses  need  to  be  made  on  turbulent 
boundary  layers  with  widely  varied  inviscid  streamline  patterns.  This 
almost  virgin  territory  for  the  experimentalist  requires  a  major  com¬ 
mitment  to  refined  experimental  technique  and  to  large-scale  facili¬ 
ties  (to  attain  the  necessary  Reynolds  number  range).  Special  atten¬ 
tion  might  be  paid  to  cases  involving  inflected  inviscid  streamlines 
and  crossover  profiles  of  secondary  velocity,  which  have  not  as  yet 
been  observed- 
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Appendtx 


NUMBERED  EQUATIONS  APPEARING  IN  TEXT 


ctu  ,  I  /'ra  +  I  ,  \  rf.„  ,01,.,  c/2,  ,  Vx  1  U  n  ,  ,  \  rU 

£  +  r-  V— 5—  +  V  ££  +  IT  (l  -  £  >  +  iKThTU5  •  2  +  x)s£ 


+  (»  +  W(i  -  e  'f ')}  -  (*)*  2-  (i  -  g'2>  -  2  (^)  *<i  -  g'l 


■  (1) 

(»_p  +  *)  £g*  +  ^  (r  +  K)(l  -£'*')+§  ^  {(•  -  f  +  >■)  gg' 

+  s(l  -  g'2)}  £'2)  +  2  x(l  -  £') 

f"'+  rx  (H1  +  »)  «' +  ir  £1  - £  ,2> +  7  ir>  {(f  +  x)  s£' 
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